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PREFACE. 

rpHE following lectures were collected from manuscripts left 
-*- by the late Professor J. C. Adams, and are now reprinted 
without change from his Collected Scientific Papers, Vol. ii., 
pp. 1—84. 

It was thought that the wide interest attaching to the 
lunar problem reached many besides the professed astronomer, 
and would justify a separate publication of this short work. 

It is known that Adams contemplated the publication of 
some such essay himself, and it must be a matter of regret to 
all that he never did so. No pains have been spared to 
present the material properly, but it is unavoidable that it 
should appear from the hands of an editor in a less perfect 
form than if the author had issued it himself 

Yet, allowing for this disadvantage, I think those best 
qualified to judge will consider this work fully worthy of 
Adams's great name. Of current elementary theories it may 
be said that they leave oflf where the difficulties of the subject 
begin, that is to say, where the various cases of slow con- 
vergence have been exposed, but not dealt with. It is perhaps 
not too much to say that these lectures carry us to the point 
where such difficulties end, in an adequate evaluation of all 
the chief constants. They leave the problem effectively solved 
and not merely stated, and shew the path clear for the for- 
mation of a detailed theory, if that is desired. 

R A. SAMPSON. 

Durham. 

8 October, 1900. 
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LECTURES ON THE LUNAR THEORY. 

[Lectures on the Lunar Theory were given by Adams 
from 1860 with few intermissions until 1889. Originally their 
aim was to illustrate geometrically the analjrtical processes and 
thereby render them more comprehensible, and they included 
some elegant theorems on the geometry of conies which have 
since become common property ; but every year several lectures 
were rewritten, and thus the whole fabric gradually changed 
into the form in which it is here presented, — the form, prac- 
tically, in which he gave them last. 

Perhaps it is superfluous to say that these Lectures stand 
upon a different footing to treatises that are intended to form 
the basis of Tables. With such, completeness is the first object 
and manner of presentation is secondary. Immense as is the 
labour of forming a treatise of this description, there exist 
several that leave little to desire in respect to fulness of 
detail. Indeed it may be suspected that their very perfection 
in the quality they profess has stifled to some degree the 
proper development of the subject, because at first sight it 
suggests that there is little left to do in the Lunar Theory, 
unless one is prepared to track down the inconsiderable errors 
that have eluded his Masters. This seems a mistake; the 
methods most suitable for the whole task adapt themselves 
comparatively ill to each detail of it, and there seems much 
that remains to be done in respect to inventing methods 
suitable for attacking separately, as far as they permit of 
separate attack, the many difficulties into which the theory 
divides at the outset, and thence perhaps approximating to 

A. L. 1 



2 LECTURES ON THE LUNAR THEORY. 

a more adequate knowledge than we now possess of the relative 
motion of Three Bodies. So far, with the notable exception 
of Dr G. W. Hill and those that have followed him, we have 
seen comparatively little eflbrt in this direction. 

This was the cardinal feature of Adams's plan, and his 
lectures shew the methods he had gradually elaborated to 
accomplish it. They separate the inequalities from one another 
as far as possible, and are content with indicating the manner 
in which these separate inequalities afterwards combine. To 
shew that, with so slight an apparatus and within so small 
a compass, the result is no mere sketch, we need but set 
side by side the coefficients of longitude found in these 
Lectures and the corresponding terms in Delaunay's TMorie. 

Adams. Delannay. 
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For those to whom the difficulties of the Lunar Theory 
are known, these numbers need no comment. 

No Manuscript exists of Lecture I. It is taken substantially 
from my own notes of 1889.] 

* With Delaunay*8 value of the Sun's ParaUax, viz. 8"*76. 



LECTURE I. 



HISTORICAL SKETCH. 



[The Lunar Theory may be said to have had its commence- 
ment with Newton. Many irregularities in the Moon's motion 
were known before his time, but it was he that first explained 
the cause of those irregularities and calculated their amounts 
from theoiy. 

Of the inequalities which are due to the action of the 
Sun, the first, — which is called the Evection, — was discovered 
by Ptolemy, who lived at Alexandria in the first half of the 
second century of our era, under the reigns of Hadrian and 
Antoninus Pius. At a very early period the relative distance 
of the Moon at different times could be told from the angle 
it subtended, and its orbit could thus be mapped out. By 
such means Ptolemy found that its form was not the same 
from month to month, and that the longer axis moved con- 
tinually though not uniformly in one direction. He represented 
this change by a motion of the centre of the ellipse, as we 
would put it, in an epicycle round the focus, obtaining thus 
a variable motion for the longer axis and a variable eccentricity. 

The representation of position by means of epicycles is 
intimately related to the modem method of developing the 
coordinates in harmonic series; thus if we have 

X = ill cos (wi^ + fli) + -4a cos (nji + 02) + •.• 
y = ilisin (rtit + ai) + il2 8in (rij* + Og) + ... 

the motion of the point (a?, y) is that on a circle of radius A^ 
with angular velocity Hj, around a centre which moves on a 

1—2 



4 LECTURES ON THE LUNAR THEORY. [LECT. 

circle of radius A^ with angular velocity Tig, and so on; and 
if, more generally, we have 

w = Ai cos (nj; + tti) + ... 

y = Bisin (nj; + ai) + ... 

we may reduce this case to *he former by rewriting 

a? = 2 (il 1 + JBi) cos (riit + oti) + g (^i - A) cos (- Wj^ - aO + • • • » 

y = 2 (Jli + A) sin(nif + ^i) + 2 i^i-^i) sin (- n^t - ai)+ ... . 

Probably we have here the reason why circular motions 
and epicycles were first employed. 

Tycho Brahe (1546 — 1601) discovered the existence of 
another inequality in the Moon's Longitude quite different 
from the Elliptic Inequality and the Evection. He found it 
bore reference to the position of the Sun with regard to the 
Moon ; so that when the Sun and the Moon were in conjunction 
or opposition or quadratures the position of the Moon was 
quite well represented by the existing theory, but from con- 
junction to the quadrature following, her position was more 
advanced than the place assigned to it, reaching a maximum 
of some 35' about half-way; and in the second quadrant it 
was just as much behind. This inequality he called the 
Variation; it was the first that Newton accounted for theo- 
retically, and if we were to suppose the Moon and Sun to 
move, except for mutual disturbance, in pure circles in the 
same plane, it is the only one that would present itself. 

The next significant step was made by Horrox (1619 — 1641) 
who represented the Evection geometrically by motion in a 
variable ellipse, and gave very approximately the law of varia- 
tion of the eccentricity and the motion of the apse. He 
supposed the focus of the orbit to move in an epicycle about 
its mean place. 

Newton's PHndpia did not profess to be and was not 
intended for a complete exposition of the Lunar Theory. It 
was fragmentary; its object was to shew that the more 



I.] HISTORICAL SKETCH. 5 

prominent irregularities admitted of explanation on his newly 
discovered theory of universal gravitation. He explained the 
Variation completely, and traced its effects in Radius Vector 
as well as in Longitude; and he also saw clearly that the 
change of eccentricity and motion of the apse that constitute 
the Evection could be explained on his principles, but he did 
not give the investigation in the Principia, even to the extent 
to which he had actually carried it. The approximations are 
more difficult in this case than in that of the Variation, and 
require to be carried further in order to furnish results of 
the same accuracy as had already been obtained by Horrox 
from observation. He was more successful in dealing with 
the motion of the node and the law of change of inclination. 
He shewed that when Sun and Node were in conjunction, 
then for nearly a month the Moon moved in a plane very 
approximately, and that the inclination of the orbit then 
reached its maximum, namely, 5° 17' about; but as the Sun 
moved away from the Node the latter also began to move, 
attaining its greatest rate when the separation was a quadrant, 
and that at this instant the inclination was 5° very nearly. 
He also assigned the law for intermediate positions. The 
fact that there was no motion when the Sun was at the 
Node, that is, in the plane of the Moon's orbit, confirmed 
his theory that these inequalities were due to the Sun's action. 

When we spoke of Newton's results as fragmentary and 
incomplete, let it not be understood that he gave only very 
rude approximations to the truth. His results are far more 
accurate than those arrived at in elementary works of the 
present day upon the subject. 

After Newton, Clairaut (1713—1765) treated the Lunar 
Theory analytically. He readily found the Variation and many 
other inequalities, but met with a difficulty in determining the 
motion of the apse. At first he made its mean motion only 
about one-half of the observed value, and supposed that this 
indicated a failure of Newton's law of the inverse square of 
the distance; but soon he recognized an error, caused by 
omission of terms which he had imagined would not affect 
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the result. When these were included the calculated amount 
was nearly doubled. 

The first Tables of the Moon which were sufficiently accu- 
rate for use in determining longitudes at sea by observation 
of Lunar Distances were those of Mayer. They obtained a 
prize oflFered by our Board of Longitude, and were published 
in 1770 by Maskelyne, the Astronomer RoyaL 

The first Theories which gave the Moon's place with an 
accuracy equal to that of observation were those of Damoiseau 
and Plana. The former was published in 1827, preceded in 
1824 by Tables; the latter was published in 1832. 

Hansen's Tables, which are those now used, were constructed 
from theory and were published in 1857 at the expense of the 
British Government.] 



LECTURE II. 

ACCELERATIONS OF THE MOON RELATIVE TO THE 
EARTH. 

When three bodies move under their mutual attraction, 
their motions are unknown to us except in the cases when 
they are approximately elliptical ; but this restriction includes 
almost all the most important cases in the Solar System. 

If one body of the system is greatly predominant and if 
the lesser bodies are not close together, the centre of gravity 
of the greater body may be taken as a common focus around 
which the others move in approximate ellipses. Or again, if 
two bodies lie close together, their relative motion may be 
approximately the same as though they were isolated, although 
the system contains a third greatly predominant body; for 
their relative motion is affected by the difference of the attrac- 
tions of the central body upon them and not by the absolute 
value of those attractions. 

The Sun and Planets are an example of the first kind; 
the Earth, Moon and Sun of the second. The Earth and 
Moon describe orbits round the Sun which are approximately 
ellipses, and the Moon might be regarded as one of the planets; 
but this point of view would not be a simple one ; the disturb- 
ing action of the Earth would be too great, though it is never 
so great as the direct attraction of the Sun, that is to say, 
never great enough to make the Moon's path convex to the 
Sun. The more convenient method is to refer the motion of 
the Moon to the Earth, and counting only the diflference of 
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the attractions of the Sun upon the Earth and upon the Moon, 
to find how this distorts the otherwise elliptical relative orbit. 
This is the method of the Lunar Theory. 

The position of the Sun must be referred to the same 
origin ; but since the Earth describes an ellipse about the Sun 
which is disturbed by the action of the Moon, if we choose as 
origin the Earth's centre, we must allow for the disturbance of 
the Sun's position by the Moon. This correction may be evaded 
by choosing as origin, not the Earth's centre, but the centre of 
gravity of the Earth and Moon, 
with respect to which the Sun 
describes a curve so closely el- 
liptical that no allowance is 
required. For, if 8y E^ M denote 
respectively the Sun, Earth, 
and Moon, and the centre 
of gravity of E and ilf, the 
accelerating forces of S are 

on E S/SE^ in ES, 

on if S/SM^ in MS; 
and these imply accelerations of of amount 

^^qr^^^]^, parallel to MS; 

now the accelerations of (S are 

E/SE^ in SE, 
M/SM' in SM; 
hence the acceleration of Q relative to /S is 
8 + E + M E „ , . p« 

S-^E + M M „ , . M-c, 

-ET-MSm'^'^^^^^^^^' 



or 



S + E + M/„ SO ^^ S0\ . ^„ 
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Let EM = r, 8Q = r\ SGM^g)] then 



E + M * E + M 



Hence 



E r 

SM^ rM E + Mr 



1 1 r -Jf r_ 

Tr: = -?; 1 ~ Tt iT- —/ " cos 0) 



and the accelerations of are 

S+E 
/I 

5 



f + Jf r EM r^. ^ ^■ au 

Now r// is approximately j^ ; neglecting the square of this 

quantity, we see that 8 moves about G in a pure ellipse. 

Consider now the accelerations of the Moon relative to the 
Earth ; subtracting the accelerations of the Earth from those of 
the Moon, we find 

let j&+if =/Lt, fif=m' ; then these become 

m'rr„ , E-Mr ( Z lb , \ 1 

parallel to OS* 
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[lect. 




In the accompanying spherical 
triangle, let Q be the centre of the 
sphere, 8M' the ecliptic, and M' 
the projection of M. 

Let Iju be the projection of 
ME on the plane of the ecliptic ; 

6 the longitude of the Moon as seen from the Earth, 

ff the longitude of the Sun as seen from G, 

8 the tangent of the Moon's latitude MM\ 

Then 

SM^to, SM'^0-0', 

r = (1 + 8^)^ u-\ cos © = cos (^ - 0") (1 + ^)-^, 

and the accelerations of M relative to E are 



1+s* 



m'(l+5«)*L E-M 1 



r*a 



1 + 



E + Mr^u 



3cos(^-^+ inif^, 



m 
? 






V 



3cos(^-^) + 
15 



E^M 1 

E + Mr'u 



(-|(l+s«) + ^cos«(^-^)) + I parallel to Gfif. 

Call these quantities U and F respectively; then if we 
resolve parallel to M'G, perpendicular to M'O in the plane of 
the ecliptic, and pei-pendicular to the plane of the ecliptic, we 
have the following quantities which we call P, T, 8; viz. : — 

r=-Fsin(^-^), 

and also 

S-P5= rscos(0-dy 

From these we find 

E'-M 1 



r, fiu^ m' ri 3 ^.^ ^. , 



E+Mr'u 



{(g-|«')co8(^-^')+^cos3(^-d')l + ...1. 
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r*M L 2 "• '^ E + Mru 

|(|-|«')sm(^-^)+^sin3(^-<9')[ + ...], 
„ „ m'sTS , 3 _,. „. , E-M 1 

{(f-|«^)co8(^-^) + ^co83(^-^)| + ...]. 

Hence with the time as independent variable we have the 
equations of motion 

d^ \dt) ' 

rdt\ dt) 

Or we may write these with as independent variable ; let 

SO that -j7 = HvJ^- 

at 

Whence H-u^{-^ + u) + H -^u'^^P; 

again, 

whence ff.«.(g + .) +F§«»| = P,-^; 
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or the equations of motion may be written 

„dH T 



ffv(g+.).p.-s-r* 



Our problem is to discuss these equations and to obtain 
from them expressions for the Moon's position at any time. 
The integration is best eflFected by observing what kinds of 
terms will disappear on substitution in the equations, and then 
assuming for the desired expressions for the coordinates a series 
of such terms multiplied by undetermined coefficients. Our 
procedure will be to discuss one by one the irregularities which 
can be isolated from one another. This will permit a survey of 
the entire field without involving needless complexity ; but if 
the Lunar Theory is to be accurate, the combinations of such 
terms with one another must also be included, and the number 
of terms employed and the labour of manipulating them becomes 
very great. 



LECTURE III. 

THE sun's coordinates IN TERMS OF THE TIME. 

To obtain the Moon's coordinates in terms of the time from 
the equations found in Lecture II., we must substitute in the 
expressions for the forces the developments of the Sun's co- 
ordinates which we now proceed to give. 

Employing as coordinates /, ff, of the last lecture, we have 
seen that the Sun's motion may be regarded as purely elliptical, 
so that 

a'_ l+e^cos(g^~iarO 

^- w' =: n'i-i!r' + 2e'sin(n'^-'Hr') + 7 e'^sin 2 (n'^ - cr')+ ... 

in which we have written for convenience v!t in place of r^t + e'. 
The quantities that enter the equations are 

Making the substitutions we find without difficulty 



(^)* 



(? 



" ^'= 1 + 1 e'« + 3e' cos (n't - w' + 1 e'» cos 2 (n't - «r') + 
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\r / sm ^ ^ "^ sin ^ ^ 



+ ^e'»'^^{(0-n't) + 2(n't-u')} 
+ 



Vr / sin ^ ' ^ ' sin ^ ^ 



+ 5<f ^1^ {3 (0 - n't) - (n't - v^')} 

-ef'^^\S{0-n't) + (n't-^')} 

+ ¥«"r{3(^-«'0-2(«'<-«r')} 



8 sin 

L ,, cos 
8 sin 



+ i^'"r{3(^-n'0 + 2(n'^-^0} 



+ 

These quantities are to be substituted where they occur in 
the expressions for the forces found in Lecture II. 

Let us now make a few general remarks upon the result of 
the substitution. 



III.] THE SUN*S COORDINATES IN TERMS OF THE TIME. 1^ 

It will be observed that the disturbing forces all involve the 
coefficient mV"'. It is very important to notice that the Sun's 
parallax is not required for the evaluation of this quantity. By 
Kepler's Third Law it is derivable from observations of the 
Sun's mean motion alone. Other terms however, namely those 
with the coefficient mja'hb, involve the Sun's parallax directly ; 
and that constant may be obtained by comparing the observed 
with the theoretical values of the coefficients of those in- 
equalities, with an accuracy probably greater than that of any 
other method. 

The mean disturbing force is radial, and is equal to 

or the mean effect of the Sun's disturbance is to diminish the 
Moon's gravity towards the Earth ; and to diminish it more, the 
greater is the eccentricity of the Sun's orbit. Now e' has been 
diminishing for ages; hence the Moon's gravity towards the 
Earth has been increasing, and its average time for accom- 
plishing a revolution about the Earth has been diminishing. 

This is one cause of the Secular Acceleration of the Moon's 
mean motion which Halley derived from the records of ancient 
eclipses. 

It may also be noticed that the coefficient of the chief 

5 
periodic part of the disturbing force, which involves 1 — h ^'^ 

increases as e diminishes. 

Finally let it be observed that the term with argument 

2(^-w'0 + 2(n'^-w'), 

which does not involve the Sun's Mean Longitude, is absent 

from the development oi \-~\ . 2(^ — ^'). 



LECTURE IV. 



THE VARIATION. 



The Variation is the first inequality we shall consider ; this 
is the inequality which is independent of eccentricities and 
mutual inclination in the orbits of the Sun and Moon. 

Let us first take the equations in the first form in which 
they are given in Lecture II., namely with t as independent 
variable : 

dt* \dtj 

rdt\ dt) ^' 

we omit the equation of motion in latitude, and in the ex- 
pressions for P, T we suppose « = 0; moreover it is possible 
and convenient to discuss separately the terms that involve 
the Sun's parallax; let these be omitted and we have 

l(?r /dey a Im'Sm' _.. .,. 



r 



and if 6'=0, / = a', wlja^^n^ 0'=^n't + €\ 

di'+rdtdi = --^n^siu2(e-n't-er, 
these are the equations to discuss. 
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Assume as a first approximation 

^ = ne + € + 62 sin {2 (n^ + e) - 2 {n't + € )} 
= n^ + e + 6a sin 2-^, say ; 

l = l[l+aaCos2i|r], 

and we shall suppose Oa, 62 so small that in the first instance we 
may neglect their squares and products. 

Substitute in the equations ; then 

4(w — 7i')* Oa cos 2-^ — [n^ + 4n (n — n') 6a cos 2^} 

+ -,{1 + 302 cos 2i|r}=in'» + ?w'« cos 2i^ 

3 

— 4 (n — n')* 62 sin 2-^ 4- 4 (w — n') wag sin 2-^ = — ^ n'* sin 2-^. 

Hence, equating the coefficients of similar terms, we have 

which gives the relation between n the Moon's mean motion, 
and - , the mean of the reciprocal of the distance ; also 

r4(n-n7 + ^la2-4w(n-n')62= |n'» (1), 

-4(n- 71762 + 4n(n-n')a2 = -|n'« (2). 

From (2) 4n (ti - n) 62 - ^^a^ = \ -^, . 

Add to (1), and substitute for /x/a'; 

3 ,, 271 -n' 



0.= 5» 



3n' - 8nn + „- n * 



A. L. 
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_ 3 ,j n(2n-n^) 1 3 n'^ 

"2"' ' (n-ny \ ^ ^ ,11 , +8(n-.nO^* 

Calline: — = m, we have 



a2 = -m* 



2-m 



2'""l-m- ^ ^11 / 
3 — 8m + -^ m* 

,3, 2-m 1 3m« 

&2 = oW^7l — -xi- n — + 



2"* •(l-w)^- o ^11 / 8(l-m)«* 
^ ^^3 — 8m + -^m^ ^ ^ 



71 m 



or, calling , , or -z = t/i,, we have 

^ w-n" 1 -m 



3 , 2 + mi 



3-2mi + 2mi« 

6.^|m.^ (^-*-^)(^;•^) -^|m.^ 

3-27^1+^7^1^ 

These are convenient expressions, and, as it happens, very 
approximate. If we wish to develope in ascending powers of 
m or mi, it appears that the latter development will be the 
more convergent. 

We find by observation — = '07480, very nearly. 

Hence 03= 00717,95, 

62= 01021,2 = 2106"-4. 

Hence the ratio of the greatest and least distances will be 

1-00717,95 : 099282,05, 

and the greatest angular deviation from the mean longitude 
will be 

35'6''-4, 

a very close approximation to the truth. 
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Also we have found 

= w2x 100280, 
which is the relation between the actual mean motion and the 
actual mean distance (or rather mean reciprocal distance) of the 
Moon. 

Without the Sun's disturbing action, the relation between 
the mean distance and the mean motion, or rather between the 
radius of the orbit supposed circular and the uniform rate of 
angular motion along it would be 

Hence in the actual orbit, the mean motion for a given 
mean distance is smaller than it would be without disturbance ; 

Or, for a given mean motion, the mean distance is smaller 
than it would be without disturbance. 

In fact, the relation between the mean distance and the 
mean motion is the same as it would be if the sum of the 
masses of the Earth and Moon were diminished in the ratio of 
100280 to 1. 



2—2 



LECTURE V. 
THE VARIATION (continued). 

We will now proceed to substitute in the differential 
equations the values of 1/r and which we have obtained, 
retaining terms of the order of the squares and products of 
Oa, 6a aiid m' or mi*. 

The values to be thus substituted are 
1 = 1(1+ a, cos 2^), 



where 



0: 


= rU + e + 


69 sin 2-^, 


^.. 


= nt+e- 


(n't+e'). 


«,= 


3 , 


2 + 7»i 


^ 3- 


-2mi + 2TO,' 



Hence 



3 

6a = (l+^)aa + ftWi'. 



• = a 1 -aaCos2'^ + ^aa'(l + cos4^) , 



^ = 4a (n — w')» [aa cos 2-^ - 2aa' cos 4-^], 
r 



dv ri 3 ~ 

^ = 4 (n - w')M 2 ^2^ + tta cos 2*^ — g aa® cos 4'«|r 
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again, 

— = n + 2 (n - n') 63 cos 2-^, 

^ J =71^ + 4^1 (n - n') 62 cos 2^ + 2 (n - n'f b^ [1 + cos 4i^], 

^ " i"{^ "^ ^"^ '''^^ ^''^ "^ I ^' (^ "^ "^^^ ^"^^l • 
Also, 



1 dr r 1 1 

-^ = 2(n-n') a2sin2^-2a2'sm4^ , 
^ ^ = 2(n— n') noasin 2'«|r + - (n - n') aA - 5 ^^a^ sin 4-^ , 



r 

\drde 
r 



— — = — 4 (n — n')* 62 sin 2-^. 

And cos 2 (^ — n't — e') = cos 2^ — 63 (1 — cos 4-^), 

sin 2 (^ — w'^ — e') = sin 2-^ + 63 sin 4-^. 

Substitute these in the differential equations, and we get, 
on transposing all the terms to the left-hand sides from the first 
equation 

4 (n — n'y ^ Oa^ + 02 cos 2-^ - ^ Oa^cos 4-^ 

-[nH2(n-n')»6aH4n(w-n')62Cos2'i|r+2(n--w')2622cos4'i|r] 

+ ^31+2^' + ^^ cos 2i^ + 2 Oa^ cos 4>|r I 

1 3 

— ^ n'2 — 5 n'^ [- 62 + cos 2'i|r + b^ cos 4-^], 

and from the second equation 

— 4 (n - ny 62 sin 2'i|r 

+ 4 (n — 71') waa sin 2'i|r + j(n - n') 0363 - 5 noA sin 4-^ 
+ ^ w'* [sin 2-^ + 62 sin 4-^]. 
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The coefficient of cos 2-^ in the first of these expressions, 
and that of sin 2-^ in the second, are respectively 

4 (n - n')' Oa - 4 yi (n - n') 6a + 3 ^ Oa - ? w ^ 

3 

and - 4 (n ~ ny b^ + 47i (ri - n') a^ + ^ n\ 

and these are evidently^ reduced to zero by giving Oa, h the 
values previously found, if we substitute for /t/a' the approxi- 
mate value n^ + a ^^''« To find the more correct value of fi/a\ 
equate to zero the constant term in the first expression ; 

2(n-w')»aa«-n,2-2(n-w')2V+^3 (l + |a.») -| w'»+|w'2 62=0, 
that is 

^3(l + |aa«) = n« + iri'«-2(ri-ri0'aa^ + 2(n-n0*6a'-|w''63 

Hence we see that fi/a^ differs from w* + ^ w'' only in terms 

of the fourth order, if we consider ttij a quantity of the first 
order and consequently Og, b^ quantities of the second order. 
Hence also by taking 

§ = n' + in'« 

in the multiplier of Og, when we equate to zero the coefficient 
of cos 2-^, we only neglected a quantity of the sixth order in mj, 
and the error in the resulting values of Og, 63 is of that order. 

We see that the substitution just made in our equations 
leaves outstanding terms of the fourth order in cos 4*^ and 
sin 4^. In order to get rid of these we must add terms of this 
form to the assumed values of 1/r and 0, respectively. Suppose 
that 

- = - [1 + o^ cos 2'«|r 4- a4 cos 4'^f'], 
0ssnt + e + b^sin 2yjr 4- 64 sin 4-^, 
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where, as we shall find, a^ and 64 are small quantities of the 
fourth order. 

It may be readily seen that the additional terms introduced 
are the following : — 

1 ePr 
in --jM 16(n — n')*a4COs4^, 

/d0 



— ( -1- j — 8n (w — n') 64 cos 4-^, 



—^ 4 3a4 cos 4-^^, 

^ -16(w-w')»648in4^, 

rd^d^ 8n(n-n)a4Sin4t, 

and also that we may neglect the terms added in the ex- 
pressions for 

|n»cos2(^-n'^-e'), | n'« sin 2 (^ - w'« - e'). 

If we write the terms thus produced along with the several 
terms left outstanding in our equations and then equate the 
whole to zero, we have 

16(n-n')'a4-8n(/i-w')64+43«4- 6 (w-n')' 02'- 2(11-71762' 

CL 

-16(n-n'>'64+8n(n-n')a4+4(»-n')*aA-2n(n-n')aj" 

from which we must determine a^, and 64. 

Put -, = ^' + hw'S 

a' 2 

g 
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and divide both equations by (n — ny\ we get 

3 / 3 \"1 27 

8(l+mi)a4-1664 + [4(l+^)-2(l+mi)]a«^ 

Simplify and multiply the last equation by ^ (1 -f- Wi), 

(19 + 6mi + 1 mi«)a4- 8 (l + mi)64-(Y+ mi -iw^^^Oa^ 

27 
-3(l+mi)wi2aa-g2^' = ^' 

(4 + 8mi + 4wi2)a4 -8 (1 + m,)6, + (1 + 27a^-\-'m^^)ai 

+ 2(l+m0(l +2mi)mi2a2 + ^(l + mi)mi*=0. 

Subtract the latter from the former and 64 will be elimi- 
nated ; we get 

(15 - 2mi + 2 rrh^a^ " ( Y + Sr^ii + | mi^^ a^* 

- ^(1 + mi) (3 +2^) ^ii'«2 - 32 (* + ^)^i' = 0, 
which gives a^ ; and this being known 64 is found from 

64 = 2 (1 + ^1)^4+ 3 (1 + ^)aa' +fg(l + 2 ^) ^^'^ "^ 256 ^'* 
Taking m = '07480 as in Lecture IV, we find 
a^ = 00004,580, 
64= 00004,237 = 8"-740. 



LECTURE VI. 

THE VARIATION [continued). 

Let us consider the problem of the Variation over again, 
taking now 6 as independent variable. 

The equations of motion are given in Lecture II : — 

^ __P T_rfw 

where _=;,_-_--_ cos 2(^-d), 

-,= -|^sin2(d-n 

so that the second equation may be written 



H' d0 



U-3n'*Qj sin 2(5-^). 



Our aim is to express t and u in terms of 6 and constant 
quantities. Now since the orbit of the Moon does not differ 
widely from a circle we may write the difference of n< + e from 
d, and the difference of au from unity as series of small periodic 
terms depending upon 0, Inspecting the form of the equations, 
it is evident that these periodic terms are of argument 2 {0 — 0') 
and its multiples ; that is 

71^ + e = ^ + periodic terms of argument 2 (^ — 0'), &c.; 

but n't + e=ff) 
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therefore 

5-^' = (l-m)^-)8 

+ periodic terms of argument 2 (1 — m) ^ — 20, &c., 
where we have written 

fi = € —me; 

this constant fi is associated with (1 — m) ^ wherever the latter 
occurs ; for brevity in writing, we shall omit it. 

We may then assume as a first approximation 

aw = 1 + tta cos (2 — 2m) 0, 

nt + €=0 + b^&iu(2- 2m) ; 
whence 

2 (^ - ^) = (2 - 2m) 5 - 2mb^ sin (2 - 2m) 0, 

cos 2 (^ - ^) = Tnia + cos (2 - 2m) - mb^ cos (4 - 4m) 0, 

sin 2 (5 - ^) = sin (2 - 2m) 5 - mftg sin (4 - 4m) 0, 

tit 
n^ = 1 + (2 - 2m) 6a cos (2 - 2m)0. 

Substitute in the right-hand member of the second equa- 
tion : — 

^^^^-^ = -3m^[sin(2-2m)^ + (2-3m)62sin(4-4m)e]. 

Therefore 

which we may write 

log^ I -J- j = 2A2COs(2 - 2m) + 2A4Cos(4- 4m) 0, 

where A is an arbitrary constant of integration, Aa is a known 
quantity, and A4 involves h^. If we take as a second approxi- 
mation 

aw = 1 + Oa cos (2 — 2m) ^ + a4 cos (4 — 4i/i) 0, 

w« + e = ^ + 62 sin (2 - 2m) ^ + 64 sin (4 - 4m) 0, 

the above value of loge{H^/h^) will not require modification 
and will supply equations of condition for determining the 
coefficients Oa, b^, cl^, 64. 
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rp, dt ^ n _ na^ h 1 

so that 

but 
loge (n ^) = -(1 -m)2 62» + (2- 2m)62Cos(2- 2m)^ 

+ [(4 - 4m) 64 - (1 - rnf 63'] cos (4 - 4m) 6, 

loge ai^ = — ^ + tta COS (2 — 2m) ^ + (^4"" "X ) ^^® (* "" *'^*) ^• 

Hence we find 

-(l-m)^6,^ = log,('^') + la,^ 

(2-2m)62 = -Aa-2a2, 

(4-4m)64-(l-m)»6a» = -A4-2a4 + |a^ 

The remaining equations of condition that we require are 
obtained from the first equation of motion ; this may be 
written 

Now 

at* = 1 + tta cos (2 — 2m) d-\- a^, cos (4 — 4m) 6, 
whence 

— ^W^ = - (2 - 2m) Oa sin (2 - 2m) tf - (4 - 4m) a^ sin (4 - 4m) 6, 

^^^ = ~(2-2m)2aaCos(2-2m)5-(4-4m)2a4Cos(4-4m)^, 

and 

(n'^Y =m3[l + (4-4m)62Cos(2-2m)5], 

(^' ^Tsin ^0-6')= m" [sin (2- 2m)5+(2 - Sm)^^ sin(4-4m)5J, 
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(n'|/cos2(^-^) 

= m^ [(2 - m) 63 + cos (2 - 2m) (9 + (2 - 3m) b^ cos (4 - 4m) 5], 

4-2 = ^ [1 + ^' " 2A2 cos (2 - 2m) 5 + (^2^ - 2A4) cos (4 - 4m) 5]. 

Substitute these in the equation above, and equate the 
coefficients of corresponding terms, 

l + |m« + |mn2-m)62+|^'aa+|m'(2-2m)aa = ^(l+A,0 

-(2-2m)X + (l + |m2V+|m« + (2-2m)m26a = ^(-2A2) 

(1 \ 3 3 

-|mn2-2m)a, = ^(A,«-2A,). 

If we neglect at first terms of the fourth order, we find from 
the first of these equations 

MCt - 1 „ 
1^ = 1 + 2"^^- . 
From the earlier set of equations we have 
(2-2m)62 = -Aa-2a2; 
substitute this in the second equation above. We get 

r-(2-2m)2 + l+|m«-2m^]a,-m%, + |m«==(l + 2m«)(--2A,), 

or 

(3-8».+^„.)«..|„.+a.=|„.+|j-^=|„.f---|, 

SO that 

3 ,2-m 1 

02 = 0^ 



and 



2 l-m^ Q.II2 
3 — 8m + -5- m^ 

h 1 1 , 

1 — m 2 — 2m 



1 3 m' 



2 



tta- 



1-m^ 8(1 -m)2 
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These are numerically equal to the quantities denoted by 
the same symbols in Lecture IV, but b^ bears the contrary sign. 
We further find 

31 



(l5 - 32m H-ymAa^ 



= I j^ {(1 + 3m - 2m») a, + (8 - 15m + Gm") b,], 

** = ~ 2 = 2^» "' - 8 '^^^ - 64 (T^ f "^ + (^ - ^^ ^^5 ' 
or reduced to numbers 

a4 = - 00002,210, 

64= -00005,414 = ira7. 

Finally let us exhibit the relation between the constants 
employed in this investigation and those of Lectures IV, V ; to 
distinguish them, attach accents to the latter, so that 

= nt + € + 6/ sin 2'*|r + 6/ sin iyft, 

— = 1 + Oa' cos 2-^ + a/ cos 4-^, 

and, omitting the constant fi as before, 

(1 - m) ^ = ^ + (1 - m) W sin 2^ + ( 1 - m) 6/ sin 4^. 

Then 

2^ = (2- 2m)^- (2 - 2m)6/sin(2 - 2m) 0, 

sin2^ = sin(2 - 2m) ^ -(1 - m)W sin (4 - 4m) 0, 

cos2i^ = (l-m)6a' + cos(2-2m)^-(l-m)6/cos(4-4m)5. 

Substitute in the equation for ; we find 

n^+.e = 5-6/8in(2-2m)^-[6/-(l-m)63Tsm(4-4m)^, 

and similarly 

a'u = 1 + (1 - m) OaV + a^' cos (2 - 2m) 

+ [a/ - (1 - m) Oa'taT cos (4 - 4m) 0. 

We observe that a' diflTers from a by quantities of the fourth 
order. 



LECTURE VII. 

CORRECTION OF APPROXIMATE SOLUTIONS. 

We may simplify the equations we have been dealing with, 
by a proper choice of units. Let the unit of distance be the 
radius of the circular orbit which the Moon, if undisturbed, 
would describe about the Earth in its actual periodic time; 
then 

[I = n\ 

Also choose the unit of time so that 
n — n' = 1, 

so that, if we take as the result of observation of the mean 
motions of the Sun and Moon, 

7i':n= 07480,13, 
we get w' = -08084,9 = mi, 

where mi is the quantity so called in Lecture IV; and 

/A =116823,4. 
We shall frequently'adopt these simplifications in what follows. 

Now let I = loga (r/a), 

so that 

rdt^dV rdt'^dt^'^Kdt)' r^' o?^ ' 
and the equations discussed in Lecture IV become 

where (o = 6 — ff. 
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Now these equations are defective, for they have been 
formed by omitting certain terms from the complete equations 
as given in Lecture II. Hence, calling lo, Oq the values of I, 0, 
which we have proved in Lecture IV to be solutions of the 
above equations, if we substitute k, 0o in the complete equa- 
tions of Lecture II, residuals are left, say X and Y respectively. 
And if I, be solutions of the complete equations, and if we 
write 

= 0,-\^S0, 

where 81, S0 are small quantities whose squares and products 
may be neglected in the first instance, we obtain the following 
equations for determining 81, 80, the corrections to approximate 
solutions lo, 00 already found : — 

Now let us write 

§ = !+„' + ., ^, = ^^-. = + .. 

where c = (1 + ny + ^ n\ 

The quantity v consists wholly of periodic terms of the form 
cos 2iy^ multiplied by small coeflScients ; w contains, besides 
periodic terms, a small constant term, which however might be 
removed if we were to choose c as the constant part of fi/vo^ in 
place of according to the definition above. 

Let 8% 8'0 be quantities defined by the equations 

then 8% 8f0 are approximations to the complete corrections 81, 
80, which if substituted in the equations that give those cor- 
rections will leave residuals, say X' and F', where 
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Z' = 2 g« ^ - 2t> ^ - 3wn + 3«'» sin 2<oS'0, 

F = 2g»^|^+2«^^ +3n" cos 20,8'^. 

We see that their value is known when 87, 8'0 are determined. 
Now B% S'0 may be determined as follows. 
Let X = po + 'ZpiCosiyp', . Y=^qi sin iyfr, 

where i takes all positive integral values ; and assume 
87 = Oo 4- Stti cos i^, 8'0 = 26i sin i-^. 
Then substituting and equating coefficients, the constant 
term gives 

Po — ScUq = 0, 

and the terms in iyfr give 

Pi - tX- - 2 (1 + nO i6i - 3cai = 0, 
qi-i%''2.(l+n')iai =0; 

the second of these may be written 



4* 

subtract from the first and we have 



% 



1i^ 

\A>x I If ^K 1 J. 



p<-2(l+«')K*-ai[t'-4.(l + «')« + 3c] = 0, 



i),-2(l+n')|-' 
or ai = s — ; 

and 6, = _2(l + n')^ + §. 

We see that a*, hi will be of the same order of small quantities 
as piy Qiy in general. And therefore the coefficients of the terms 
of X\ T* will be of order higher than those of X, F. Proceed 
then to determine further corrections 8"Z, h"d satisfying the 
equations 

F-.*^ + 2<l + »0^ -0; 
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then if S7 + S'7, S'0 + S''0 are substituted in the complete 
equations for SI, B0 the residuals become 

X" = 2f^^ -2v^^ -SwS"l + ^'*8m2c,^'0. 

F" = 2§^ + 2v^4? +3n'»cos2«S"^, 

at at at 

expressions which, if developed in series of cosines and sines of 
multiples of *^, will have coefficients of higher order than the 
corresponding coefficients in X\ Y\ The like process may be 
repeated until the residuals become insensible; we then have 
sensibly connect values of SI, S0, giving 

l = lo + Sl = k'^S'l + Sl + , 

0^0o + S0=^0o + S'0+S''0 + 

We may now take into account squares and products of the 
small quantities SI, S0 by treating lo + Sly 0o + S0 as given 
approximate solutions just as we have here treated l^, 0^; 
substitute them in the complete equations of motion, and 
determine the residuals X, Y which they leave. These resi- 
duals will form the basis of a second approximation, and the 
operation may be repeated until no further correction is 
necessary. It is to be observed that if SI, S0 depend upon 
some such constant as the eccentricity of the Earth's orbit 
around the Sun, or the parallax of the Sun, then successive 
approximations yield correctly and separately the terms which 
depend upon the first, second, powers of that constant. 



A. L. 



LECTURE VIII. 

THE PARALLACTIC INEQUALITY. 

We shall now apply the method of the last lecture to find 
the terms in the Moon's coordinates which depend upon the 
parallax of the Sun. 

The values of I, 6 found in Lecture IV are 
/o = \oge (r/a) = - Oa cos 2^, 
^0 = wf + 6 + ftasin 2^, 
and these satisfy the equations of motion in which the terms 
involving the Sun's parallax are omitted. Hence the residuals 
they leave from the complete equations are 

Z = -X/i'«^|^cos(^o-^) + ^cos3(^o-^')}. 

Y= Xn'«^||sin(^o-^) + ^sin3(^o-^')}, 

, ^ E-Ma 

where X= ^ — r^-;. 

Now from above 

^,_^' = ^ + 6asin2i^; 
hence we have 

sin(^o-^)= sin^ + ^ftaCsin^ + sinS^), 



cos (^0 - 0') = cos -^ - - 62 (cos -^ - cos 3-^), 
sin 3 (^0 - 0') = sin 3V^ + 2 62 (sin -^ + sin 5>^), 
cos 3 (^0 - ^) = cos 3^ - ^ 62 (cos '^ - cos 5-^); 
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and 

^ |?cos(^,-^)+^co83(^. - ^)} = (|-^ 6, - |a,) co8^ 
^ /15 , 9 . 9 \ o I . /*5 , 15 \ , , 

/16 3 . 3 \ . „ , /45 , 15 \ . , , 

Assume 

— 8\ = Xtti cos yfr + A^s cos 3*^, 
S^ = XAi sin ^ + AAs sin 3*^, 
neglecting for the present the terms in 5-^. In the present 
case it happens that it is more advantageous to substitute 
these expressions directly in the complete equations for SI, S0 
given in the last lecture than to follow exactly the process for 
finding them by successive approximation. Omitting the 
factor X, we get 

Oi cos '^ + dos cos 3-^ + 402 sin 2^ [aj sin sft + 3a, sin 3^] 

+3 -^[aiCOS'^+asCosSV^] +3 -3[3aaCos 2-^] [aiCos'^+ OjCOsS-^] 

-[2(1+ n') + #2 COS 2yjr] [bi cos yjr + Sb^ cos 3i^] 
-♦- 3n'" [sin 2-^ + 62 sin 4-^] [bi sin yjr + b^ sin 3-^] 

,,/9 27 , 3 \ , , /15 9 , 9 \ ^, 

,/45, 15 \ _ ^ 
-nM jg 62 - jg aajcos6Vr = 0, 

— 61 sin -^ — 96s sin 3-^ + 4aa sin 2-^ [61 cos -^ + 36, cos 3^] 

+ 3n'" [- 6, + cos 2-^ + 62 cos 4-^] [61 sin ^ + 6, sin 3-^] 

+[2 (1 + n) + 46a cos 2^] [oj sin >^ + 3as sin 3-^] 

. ,,/3 21 , 3 \ . . , /15 3 , 3 \ . ^, 
+n«(3-g-62-ja2JsinVr + n»(-g- + ^62-^a2Jsm3Vr 

/45 15 \ 

3—2 
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If we equate to zero the coefficients of cos yjr and cos 3-^ in 
the first, and those of sin yfr and sin 3-^ in the second, we obtain 
the following equations for aj, 6i, Oj, 63 ; the terms in 5-^ remain 
outstanding, and the effect of ag, 65 in modifying the other co- 
efficient is neglected. 

a,[l + 2a, + 3^,(l+|a,)]-6,[2(l + n') + 26,-|n'^] 
+a,[6a, + |^, J-6, r66,-|n"(l+6,)l = |n'«(3-96,-4a,), 

a, [2 (1 + »') - 26J - 6, [1+ 1 «'» - 20, + 3n'»6,l 
+a, [66J -6, r6a,-|n'»+|n'«6al = -|n'»(l -76,-20,), 

+a,[9 + 3^,] -6,[6(l + n')] = |«"(« + |*«-1«»). 

+ a,[6(l+n')] -6,[9 + 3w'»6J |'*''(^+i*»-i°»)' 

If we require the formal values of Ui, 6„ Oj, 6„ we must 
substitute for a,, 62, /*/a' the expressions we have found for 
them, and it will then be best to develope the coefficients in 
ascending powers of n'. But it is difficult to obtain by this 
process such good numerical results as we can get by sub- 
stituting the numerical values of a,, 6j, /*/a' immediately in 
the equations above. If we do so we get the equations 

4o6672a,-2-172326i+ 080930,- 061376,= |ro''x 2-87937, 
214128a,-0-995646,+ 061270,- •033386,=-|w''x 0-91416, 

O 

-02349ai- -030126i+12-61461a,- 6-485086,= |n'»x 600455, 

O 

•02042aj+ -024066,+ 6485080,- 9000206,=-|w'»x 500152. 

O 
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We notice that the first equation is not very different from 
the second doubled ; it is this fact that makes successive 
approximation a disadvantageous method and renders it ad- 
visable to include small quantities from the beginning. 

Eliminate a,, 6, in succession from the third and fourth 
equations, thus : — 

Multiply the third equation by 
9-00020-r[12%514olx 900020 6-48o08x6-48508]= 0-127523, 
and the fourth by 

-6-48508-r-[12-51451x9-00020-6-48508x6 48508]=-0-09]887, 
and add ; 6, will be eliminated. 

Again multiply the third equation by '091887 and the 
fourth by —•177317 and add; Og will be eliminated. Hence 
we find 

•001119ai--0060526i4-a8=in''* x 1*09776, 

- -0014630, - 0070346, + 6, = | n'» x 1-34669. 

Multiply these by —'08093 and '05137 respectively and add 
to the first equation : — 

4'o6672a,-2'l72326,+ 08093a,-051376,=|«'«x 2'87937, 
-0'00009a,+0'000496,-'08093a, =|n'»x -0'08884, 

O 

-0000080,-0000366, +0ol376,=|»'»x 0'06918; 

O 

hence 

4-5666501 - 2-l72196i = I n « x 285971. 

o 

Eliminate a^^ 63 in a similar manner from the second 
equation ; 

2-14128ai-0-995466a+-06127as--033386,=-|n'2x 091416, 

o 

-000007a,+0-000376,-'06127o, =-|n"x 006726, 

O 

-0'00005a,-0'000236, + '033386,= -|?i'» x -004495 ; 

O 
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hence 

214116ai-0-995506i=-|n''' x 093647. 
From these equations we find 

a, = -|n'«x46-4814 =- -11392,8, 

o 

6, = _ I n'» X 99-0336 = - -24273,4, 
0,== |«"x -56042= -00134,9, 
6,= |»'»x -58209= 00142.7. 

O 



LECTURE IX. 

THE PARALLACTIC INEQUALITY (continued). 

Let us now consider the terms in 5-^ which have been left 
outstanding. 

Include additional terms Tat^ cos 6*^, Xb^ sin 5^ in — Bl, BO, 
and equate to zero the coefficients of cos 5^, sin 5-^ in the 
differential equations that give SZ, BO. 

We have 

2oa, + ^a, - 10 (1 4- n')b, - 6a^, - 66A + s ^ ci^z 



-256, +10(l+n')a»+6as6, + 6a,6s 

3 ,„, . „/45, 15 



is + 60362 



+ 2 



In these equations substitute 



g = l + 2.'+|.'. 



Then 



(28 + 6w' + I «') a. - 10 (1 + n') 6, = (I - 9»' - ^ w'») OaO, 
-256, +-10(l+»')a.=-6a,6,-66^- |n'»6,-n'>(^6,-^«,) . 
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Eliminate 65 : 
(24 - 2n' + I n'«) a. = [(I - 9»' - ^ »'•) a, + ^ (1 + n') ft.] a, 

+ [66,+ |n" + ^|(1 + n')a, + |(l + »')«'*] h 

and then 6, is given by 

From these we find 

o, = I n'» X 00595,3 = -00001 4591, 

6. = |«'» X 00710,3 = 00001,7410. 

O 

These numbers being so small, we see that we may safely 
ignore, as we have done, their effect in modifying the earlier 
coefficients. 

To find the effect of these coefficients upon the Moon's 

coordinates we must multiply by the factor \ = p, — jj^ . — . 

We shall take in accordance with the results given in 
Monthly Notices, Vol. 13, p. 177, and Appendix to the Nautical 
Almanac, 1856, 

1 = 81-5. 

Constant of Moon's Parallax = 3422"-325. 

Also .we shall take in the first place, the Sun's Mean 
Parallax to be 8-':8, and in the next place 8"*9, and we will 
find the corresponding values of the coefficients of the Paral- 
lactic Inequalities. 
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We find 



8"-8 



8"^9 



\= 


•00250,9 




\ = 


•00253,76 


}Mi = - 


-•00028,585 




Xa. = - 


-•00028,910 


X6. = - 


-00060,903 =- 


- 126"^62 


X6i = - 


- 00061,696 
= -127"05 


Xa8 = 


•00000,3385 




XOjSs 


•00000,3423 


XA.= 


•00000.3680 = 


0"7384 


X6.= 


•00000,3620 

0"-7468 


\a.= 


•00000,00366 




Xa.= 


•00000,00370 


X6.= 


•00000,00437 = 


0" 00901 


XA.= 


•00000,00442 

0"00911. 



These results are very fairly accurate ; but in order to get 
good values for «!, 6a, we were obliged to discuss Oi, 6i, (hyh 
simultaneously. Let us consider the peculiarity of the equa- 
tions from which this difficulty arose. 

Following the method of approximation of Lecture VII, if 
we neglect at first the products of 8/, BO, dBl/cU, dSO/dt with the 
small quantities 02* ^j, n\ the equations become 



d^ 
d'B0 






dt 



W'^^^'-dt 



+ F = 0. 



Now suppose the following is a set of terms that appear 
in X Pi CO8 (it + y)y in Y qi sin (it + y), 

Bl Oi cos (it + 7), B0 bi sin (it + 7) ; 

then as in Lecture VII, we find 

Pi-2^qi 



ai = * 



i^-n'-h 



3 
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Therefore if % differs little from w, the divisor in Oi will be 
small, and a small error or omision in the numerator of Oi will 
appear magnified in the values of both Oi and ft*-. In the case 
of the first term of the Parallactic Inequality, 

i = ?i — n\ 
f»- n«4- 1 w'2 = - 2nn' + |n'* ; 
and if we take 

which differ from the correct values by quantities of the fourth 
order, then 

and the formulae give 

-? w'(o/i-3n') 
^*-4(n-V)(4n-5n'r 
, 2n 3 n'« 



Now if we develope these expressions in ascending powers of m, 
i.e. n'/uy the first terms are 

15 , 15 

and these are the only terms which the formulae derived from 
our method of approximation will give correctly. 



LECTUEE X. 



THE ANNUAL EQUATION. 

Let us next take into account the effect of the first power 
of the eccentricity of the Earth's orbit. We shall find that it 
produces an inequality in the Moon's coordinates, the chief part 
of which has a period of one year, and is therefore called the 
Annual Equation. 

In the formulae of Lecture VII, let the known approximate 
solutions Iq, 0Qy include the Variation only ; then the equations 
for the corrections Zl, SO are 

X-h^/-h4a2sin2>f^^^-3£3(H-3a2COs2Vr)8i 

- 2 [(1 + n') + 262 cos 2f ] ~ + 3n'«(sin 2^+ b^ sin 4^^)80 = 0, 
7+ ^ + 4a^ sin 2i/r ^ + 3n'n- i>2 + cos 2>^ + 62 cos 4Vr) S^ 
+ 2 [(1 + n') + 26, cos 2ylr] ^ = 0, 

where a^, 62, /tt/a* are known quantities whose values are given 
in Lectures IV, V. 

Refer now to Lecture III, and we find that the terms that 
are left outstanding when the terms of the Variation are sub- 
stituted, and the parallactic terms omitted are the following : — 



44 LECTURES ON THE LUNAR THEORY. [LECT. 

Z = - 1 wV COS in't - w') - ^ ^''«' COS {2(^ - n't) - {n't - vr')] 

+ I nV cos {2 (^ - n't) + (n'« - m)] 
F= + ^ nV sin {2 (^ - n't) - (n'^ - vr')] 

- I /iV sin {2 (^ - w'O + (n'^ - ^% 

Write a for n't — ur' ; then 
cos {2 {0 - n'O - «} = cos (2-1^ - a) - (262 sin 2-^) sin {2-^ - a) 
= - 62 cos a + cos (2^ — a) + 6a cos (4^ — a), 
sin )2 (^ — n'f) — a} = 4- 62 sin a + sin (2-^ - a) + 62 sin (4^ — a), 
cos {2 (^ - n't) + a} = — 62 cos a 4- cos (2-^ + a) + 62 cos (4^ + a), 
sin {2 {0 - w'f ) + a} = - 62 sin a + sin (2i|r -f- a) + 62 sin (4^ 4- a). 

Hence 

Z = -|w'2(l-362)e'cosa-^nVcos(2>f^-a) 

3 21 3 

+ 7 nVcos (2i^ + «) - x w''62e'cos (4i^-a)+ 7 w'^ftge' cos (4i^-»-a), 

F = 6w'%2e' sin a + ^ nV sin (2Vr - a) - 1 w V sin (2^ + a) 

21 3 

4- -J- ^'*2e' sin (4^ - ^) - 4 '^'^^^' sin (4^ 4- a). 

For our present purpose we shall ignore the small terms in 
4*^ — a and 4*^ 4- a which are of the sixth order. 

Assume 

— Si = a^e' cos a 4- a^fi' cos (2^ — a) 4- (he' cos (2'i/r 4- a), 

S^ = 656' sin a 4- ftse' sin (2-^ — a) 4- 676 sin (2^- 4- a). 

Now the terms which arise in the left-hand members of the 



X.] THE ANNUAL EQUATION. 45 

equations owing to terms —Op cos pt in Bl, and bp sin pt in BO, 
will be 

p^Op cos pt + 2aapap [cos (pt — 2-^) — cos {pt + 2-^)] 

+ -^ dpi cos pt + '^(hcosi pt — 2'i/r) 4- 2 «« cos (^^ + 2'i/r) 

- 2(1 4-w')i>ftpCosp«- 262i?6ii[cos(^^- 2^)+cos(^^+2'^)] 

+ 1 n'% [cos (;)^ - 2^) - cos (pt + 2^)], 
and 
— p\ sin p^ + ia^pbp [— sin (p^ - 2-^) + sin (pt + 2*^)] 

Q 

+ 2 n'% [sin ( p^ - 2-^) - 263 sin pt + sin (pt 4- 2>^)] 

+ 2 (1 + n')pap sin pt + 262pap [sin(p^ — 2'^)+sin(p^+ 2^)], 
respectively, neglecting the very small quantities in 4-^. 

Hence we get the equations following : — 
Equate to zero the coefficients of cos a, sin a : 

- 2 (2 + «') ^A + 1 M'^ftr = I «'' (1 - 36,) 
- n'% - 3n'»6A + 2(1+ n') n'a, + 2 (2 - n') ajb^ 
-2(2-n')ft^.-|«''6.-2(2 + n')aA 

+ 2 (2 + n') 6,0, + 1 n''6, = - 6n'%. 
Equate the coeflScients of cos (2'^ — a), sin (2'^ — a) : — 
(2 - n')»a, + 3 ^, a, - 2 (1 + »') (2 - wO ft* + 2n'a^. + 3 ■^, . I OiA. 

-2n'hA + %n'%= ^n'\ 
- (2 - nj 6, - 3n'»6A + 2 (1 + »') (2 - n) a, + 2»'a,6, - 2n'h4i, 
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Equate the coefficients of cos (2^ + a), sin (2-^ + a) : — 

Ot Ob M 

- (2 f n'f br - 3»'*6A + 2 (1 + n') (2 + n')a, + 2n'a^, + 2n'b^ 

+ |«'»6.= In- 

In equations of this class, as a general rule we would 
determine ag, 65 approximately from the first pair, substitute 
them in the second pair and determine a„ 6^ approximately, 
and similarly 07, by from the third pair, and repeat this approxi- 
mation as often as might be necessary. But if we refer to the 
second equation, we see that 65 must be determined by means 
of a small divisor, and this puts any method of approximation 
at a disadvantage. In order to obtain readily satisfactory values 
for the new coefficients, we shall treat the six equations simul- 
taneously, substituting first the numerical values of the known 
quantities. 

We have found 

a, = 00717,95, b^ = 01021,20, /i/a» = 1-17150,3. 
Hence 
3-52105a5 - 017476365 + •065406a« - -02939365 + 06772707 

- 03269567 = I w'2 X 0969364, 

017476305 - 000673665 - 03919706 + 01775366 + 04249907 

- 02007567 = I ^'' X - 0040848, 

•O39OO905 + 00815365 + 71976606 - 41485866 

= |n'2x 3-50, 



- 00165105- 00864365 + 4-1486806- 3-683356, 



'6 



Q 
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•036687a.- •0114566. +7-84443a,- 44981 16, 

= |«'»x-0-50, 
•001651a, + •0109656. + 449811a, -4^3301 26, 

= |n''+ 0-50. 

From the second and third pairs we find 

•016186a, + 0070846, + a, = | «'» x 2-94739, 

•018678O, + •0103266. + 6, = | »'» x 4 26994, 

■011026a, - •0072036. + a, = | »'» x - 0321398, 

•011072a. - 0100156. + 6> = | «'» x - 0449338. 

Eliminate a«, 6„ a,, 6^ from the first pair. 
Hence 

3-52015O. - •1747626, = | «'' x 0909170, 

0174819a, - -0065366. = 5 n'' x 0003516. 
Hence 

a, = I «'» X - 070619,8 = - ^00692,37, 

6. = I w'» X - 19-4268 = - 19046,3. 

Now e' is a -constant found by observation; taking 
e' = 3459"'28, its value in 1850, we get 
a/ = -00011,61, 

6,e' = - -00319,4 = - 658"-9. 
and further 

a,= -03035,8, a,e'= 00050,9, 

6,= 04396,7, 6.«'= •00073,73= 162"-09, 

a, = - •00444,7, a,«' = - -00007,457, 

6, = - 00623,6, 6,«' = - 00010,46 = - 21"^57. 



LECTURE XI. 

THE EQUATION OP THE CENTRE AND THE EVECTION. 
We have seen that the equations of motion 

are satisfied very approximately by the values 

= n^ + € + 62 sin 2-^, 

where yjr==nt + €- (v!t + e'), 

and a-2, 62 are small quantities depending upon the ratio w'/n, 
and a is a quantity depending upon n in such a way that 

while w, € are arbitrary, though subject to the assumption that 
the ratio njn is small. 

This solution, then, expresses a possible case of motion; 
nevertheless it is no more than a particular case because it 
involves only two arbitrary constants, whereas the complete 
and general solution must contain four, in order that it may 
be able to satisfy any given initial conditions, that is, in order 
that the initial coordinates and their initial velocities may have 
any given values. 
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When there is no disturbance the four arbitrary constants 
are n and €, — which denote quantities similar to those expressed 
by the same symbols above, — and the two elliptic elements e 
and isr, of which e denotes the eccentricity of the orbit and 
w the longitude of the apse. 

We will now shew how to complete the solution by intro- 
ducing into log(a/r) and additional terms depending on 
quantities similar to e, isr, of which the former ife constant and 
the latter varies slowly and uniformly with t ; and for the sake 
of simplicity we will suppose at first that e is so small that its 
square and higher powers may be neglected though it is other- 
wise arbitrary in magnitude. 

Let us assume then 

log - = Oa cos 2yfr -h e COS (nt — tsr), 

^ = n^ + € + 62 sin 2-^ + 26(1 + 60) sin (n^-tsr), 

in which the elliptic terms are of the same form as in the 
undisturbed orbit, and w is supposed to be slowly variable, 
so that 

where p is supposed to be a small quantity of the order of the 
disturbing force. 

We will now substitute these assumed values in the differ- 
ential equations. We have 

-Tj: = 2 (n — n') Oa sin 2-^ + (n—p)e sin (nt — isr), 

dH 

-^ = 4 (n — ny Oa cos 2^ + (w — pfe cos (nt - w), 

— = n + 2 (w — 7i') 62 cos 2-^ -}- 2 (n — p) ( 1 + 60) ^ cos (nt — -Br); 

_ = — 4 (n — n') 62 sin 2-^ — 2 (n — ^^ (1 + 60) ^ sin (nt — w). 

Hence 
4 (n — ny a^ cos 2-^ + (n — pY e cos (nt — ot) 
+ 2 (w — w') (n''p)a^ [cos (2^ -nt^-vr) — cos (2-^ + w^ — -bt)] 
A. L. 4 
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— {n* + 4n(w -n') 62 COS 2-^ + 471 (n — p)(l +60) « COS (n^ — OT) 

+ 4(n-n')(n-p)(l+6o)62c[co8(2^-r?f+'Br)+cos(2^+n^-'Br)]} 

+ -3(1 + Soa cos 2^} {1 + Se cos (nt - -or)} 

- n'^ I2 + 2 cos 2^-3 sin 2^ [2 (1 + 60) e sin (w^ - bt)]! = 0, 
and 

- 4 (ri - w')^ 6a sin 2-^ - 2 (n -^y (1 + 60) e sin (n« - -or) 
+ 4w (n — w') Oa sin 2-^ + 2w (n — p) e sin (nt — isr) 

+ 4(n--7i')(w— p)(l+6o)««a[sin(2'^— w^+'Br)+sin(2'^+n^— -or)] 
+ 2 (n - w')(w — p) e6a [— sin (2-^ — n^ + 'Br) + sin (2-^ + rU — -bt)] 

+ nf^ L sin 2Vr + 3 cos 2-^ [2 (1 + 60) e sin (w« - "bt)]! = 0. 

It will of course be found that with the values of Oa, 6a of 
Lecture IV, the terms independent of e vanish identically. 

Equating to zero the coefficients of cos (nt — -or) in the first 
equation and sin (nt — isr) in the second, we get 

(n-.p)2~4ri(n-p)(l + 6,)+?^ = 0, 

- 2 (w -p)2(l + 60) + 2n (n-p) = 0. 

Therefore (n — p) (1 + 6©) = w, 

and (n-py=4n2--^ 

3 

=^n^ — - w'^ approximately, 

p 3 

or - = T ^'^ = i!>o> approximately. 

Now terms have been left outstanding with the arguments 
2-^ — 71^ -h OT, 2-^ + n^ — "BT. These may be removed by assuming 

log - = Oa cos 2-^ + e cos (n^ — tsr) + a^e cos (2'\^ — n^ -I- -bt) 

+ Oaa^ COS (2-^ + n^ — "or), 

^ = n^ + € + 6asin2'<|r + 2e(H-6o)8in(7i«-'Br) 

+ h^e sin (2^ — n^ -h ot) + 6aae sin (2^ -f n^ — «r). 
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Hence in place of the former equations, we get the 
following 

(„__p)._+„(„_p)(l + 6.)+^ 

+ r2(n-n')(n-2n' + p) + |g1a,aa 
+ h(n-n'){Sn-2n'-p) + ^^^a^a^ 
+ [-2(n- n') (n - 2n' +i)) 6, + 1 n'^l b^ 

+ [-2(n-n')(Sn-2n'-p)b^ + ^n'*]b^=0, 

-2(n-py(l + b,) + 2n (n-p)- 6n'»6, (1 + 6.) 

- 2 (n - nO (n - 2n'+p) b^a^ + 2 (« - «') (3n - 2n' -p) 6,0^ 

■\-[2(n-n')(n-2n'+p)a,-^n'»\ka 

+ r-2(n-n')(3»-2n'-p)a, + |n'»l6a = 0. 

2n 
Multiply the second by -— , and add to the first; this 

will eliminate 1 + &o- 

+ 2 (n - n) (n - 2n' +p) [a^^ - 62621] 
9 a 3 

+ 2 (n - n') (3n - 2n' -;>) [o^Oaa - b,b^] 

it ~~ X? 

- 4 (n - n) (3w - 2w' - p) [63022 - 02622] 

Tl *~" p 



n— p"- ^ 



4—2 
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Also the equations obtained by equating the coefficients of 
e cos (2'\^ — n^ + ot) and e sin (2-^ — n^ + w) to zero are 

4 (ri — n') (n -p) (1 + 6o) <3t2 - 2 (n — v!) (n —p) 63 

- 3n'2 (1 + 60) - (?! - 2n' +p)2 621 + 2n (n - 2ri' + i^) Oai = 0. 

2n 
Multiply the second by — _«-/ , and add to the first; 

this will eliminate 621 > and gives 

[2(n-«')(n-p) + |g-„--|^(n-«')(n-i))(l + 6.)]a, 

Lastly the equations obtained by equating the coefficients of 
e cos (2-^ + n^ — -or) and e sin (2^ + nt — m) 
to zero are 

-2(n-n0(»^-;>)+|g]a.-4(n-n')(n-p)(l+6o)62-3n'>(l+6o) 

+ r(3w-2n'-;>)« + ^^1a22-2w(3w-2n'-;))622==0, 

and 

4(n-n')(ri-;))(l + 6o)a2 + 2(n-r2')(w-;))62 + 3n'>(l + 6o) 

-(3n-27i'-;?)»622 + 2n(3n-2n'-p)aa = 0. 

2n 
Multiply the second by — ^ : _o ^_ ^^^ ^^^ *^ *''^® ^^^ > 

this will eliminate b^, and gives 
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+ \(Sn - 2n' - pf - 4?i-2 + ^^1 o^a = 0. 

These six equations are to be solved by successive approxi- 
mation ; taking the first rough values of p/n and 6o> we find 
from the last two pairs values for a^, 621* <^22> i!>22; these are 
substituted in the first pair and 3deld more approximate values 
of p/n and 6o» ^^^ so on. 

It will be noticed that this complexity is made necessary by 
the fact that Oa, 621 are found by means of a small divisor 

(n-2n' + py-4fn:'-{-^. 



LECTURE XII. 

THE EVECTION AND THE MOTION OF THE APSE. 

We proceed to the conversioD into numbers of the formulae 
of Lecture XI. 

Take n — w' = 1, 

n= 1-08084,9, 
ri'= -08084,9, 

§ = 1-17150,3, 

log aa = 7-85609, 
log 62 = 8-00911. 

First Approadmation, 

(n-p)^-^4n« + 3^ = 0, 

4^2 4-67293,7, 

- 3 ^, - 3-51450,9, 

(n-py^: 1-15842,8, 

n-p= 1-07630,3, 

p = -00454,6, 

n-2n'-\-p= -92369,7, 

3n - 2n' - ;) = 3-07630,3, 

1 + 60 = 1 00422,4 = n/{n - p). 
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Substitute in the equation for a^, 



The various terms give 

2(n- n') {n-p)a^ '01545,45 

I^^Oa -03784,83. 

- 8 J^ (n - n') (71 - ») (1 + bo) (h - '07264,10 

- 4 (w - n') (n -p)il+ bo) 6, - '04415,05 

4 — ^7 (w - n') (n - _p) b, '05144,47 

3n'2(l+6o) '01969,25 

6 ^Ct- (1 + ^o) -04608,58 

•05373,43 

(n-2n'+py -85321,6 

-4n2 + 3^, -115842,8 

a' 

- -30521,2 



031 = -17605,6. 

The equation for 621 is 

1 2n . n — n . v/^iv 

621=?= / . 021 + 4;; ^r-r-; — xi(^-p)(l + Wa2 

n — 2n 4-^ (71 — 271 +py ^ ^ 

71 -7i' ^'2 



(>.-2n'+^ )^^"-^P>^-- ^»-2»+j,)' <^ + ^')^ 
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271 



a^ -41201,8 

•03637,95 



- 2 . \~ "' ^ (» - p) i, - -02576,42 



6a = -39965,3 
The equation for 0,2 is 

-83;p|^(n-n')(n-l')(l+W«.-4(«-»')(n-i')(l+^.)&2 

-«3;^^i:^<i-^^«)=«- 

Here - 2 (w - ri') (n - jd) «« - '01545,45 

l^.o, 03784,83 

-8g^^^J^^(n-n')(»-;))(l+6.)a. --02181,13 

- 4 (n - »') (n - !>) (1 + 6.) h - -04415,05 

'Sn-2n'-p^ 

-3n'ni + 6o) -'01969,25 



- 4 ^^_ol'_,. (»* - «') (" -i>) *» - -01544,69 

6„) -'01969,25 

(l+6o) --01383,78 



3»-2w'-p 



- -09264,52 



(■^-2n'-py 9-46363,4 

-4»'+3^. -1-15842.8 
a' 

8-30520,6 
ttjg= 01114,30. 
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And 

71 -^ W 71 ^ 

+ ^3«-2n--p)' <"-^>^'-^^ 3«-2n--p)' <^ + ^'>' 
4 (3^ "^/_^), (»^-.p)(l + &.)«. -00327,988 

^3n-2"n--;>/ ^-^^'> '^^^^Q^'^^^ 

6a = -01651,37 



Second Approximation. The complete equation forn — ^ is 

(n -!>)» - 4«' + 3 ^, + 12 — 6, (1 + 6,) 

9 11 

+ 2 (n - n*) (« - 2n' + 1)) [ojOa - 6j6j,] + 5 -, Kon + aaOu] 

+ 2 »" [^n + 6»] + 2 (» - n') (3n - 2n' - p) [0,0^ - 6,6 J 
+ 4 — ^ (n - w') (m - 2»' - !>) [6,o„ - a,6a] 
-4 (7i-n')(3n- 27i' — p) [62022 — 02622] 

71 "^ 2J 

-4»» + 3^, -115842,8 

12—6,(1 + 6.) -00080,8 

2(w-n')(n-2n' + p)[a,aa-6,6a] - -00520,1 
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l^A'^<hi + (ha^] -00708,4 

|n'»[6a + 6J -00406,9 

2(«-n')(3n-2n'-p)[a,aa,-M«] - '00048,2 

4,-^(n-n')(n-2n'+p)[b^a^-a,b^] - -00397,5 

- 4>-^{n-n')(3n-2n' -p)[btaa- (hb„] - -00003,0 

" ""''■[t^i-tj -00756,1 



n—p 



(n-p)*= 1-14859,4 
n-p = 107172,5 
p= 00912,4 
p:H= -00844,2. 

Apply these numbers in the equation for 1 + 6, : — 
l+*»=^-3(i^&.(l+t.)-(-^(«-2»'+i>)[6.a.-a>6„] 



n 



n— /) 



1-00851,33 



-^(ir=^M'^ + bo) - -00017,51 



w — n 



t;;^ y (« - 2m' + p) [6,Oa - a,6a] 



•00086,56 



^^-^(3n-2re'-p)[6,a„-a,6J -00000,64 



l+6o= 100757,1 



Continuing the approximation for a^, ^n, a„, ^a tl^e various 
terms found are the following : — 
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•01638,88 
•03784,83 
•07221,53 
•04410,93 
•05097,33 
•01976,81 
•04601,14 
•06365,63 



divisor 
•86169,6 
- 1-15842,8 
- 29673,3 



aa= 18082,0 



- 01638,88 
•03784,83 

- -02182,34 

- •04410,93 
-01640,41 

- ^01975,82 

- 01390,46 
-09264,01 



divisor 
943649,4 
- 11.5842,8 
8^27706,6 



aa = -01118,03 



•42108,0 
•03598,79 

- 02540,22 

- 02292,94 
b^= -40873,6 



•00786,805 
•00328,669 
•00231,985 
•00209,403 
62, = •01556,85 



Third Approanmation. We find 



115842,8 

- 00081,40 
•00535,02 

- -00723,64 

- -00416,02 
-00048,35 
-00407,40 
-00002,97 

- -00777,57 
(n-py= 1-14838,0 



n-|) = 107162,5 

p= -00922,4 

p:n= -00853,5 



-01538,73 
-03784,83 

■ 07220,64 

■ -04410,81 
-05096,30 
-01975,95 
-04600,96 

•05366,42 aa = ^18092,9. 



divisor 

•86188,04 
- 1-16842,8 
- -29664,8 



1-00860,75 

- -00017,57 
-00087,94 
•00000,64 

- ^00167,84 
1+6.= 100763,9 



-01638,73 divisor 

•03784,83 9-43488,0 

-•02182,35 -1-15842,8 

-•04410,81 8^27645,2 

-01640,32 

- 01975,95 
-01390,60 

- 09253,93 aa!=01118,10. 



•42128,90 
•03597,92 

- 02539,43 

- 02292,60 
ba= -40894,8 



-00786,881 
-00328,671 
•00231,978 
•0 0209,430 
6,, = •01556,96 
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Fourth Approasimation. 



n-p = 1-07160,3 

p= -00924,6 

p:n= 00865,4 



1-00862,79 

- -00017,67 
-00088,04 
•00000,64 

- 0167,94 
1+6,= 1-00766,0 



1-15842,8 

- -00081,41 
•00534,22 

- -00727,13 

- -00416,24 
•00048,35 
-00407,66 
-00002,97 

- -00778,05 
in-py= 1-14833,2 

The values already found for the remaining quantities are 
sufficiently exact. 

These numbers give, taking after Hansen, 
e (1 + 6,) = -06491 
e = -06449 
2 (1 + 6,) e = -10982,0 = 22651"-9 
Oae= 00986,03 
6aC = -02228,44 = 4596"-6 
Oae = -00060,93 
6^6 =00084,85= 176"-1, 

and taking the Moon's mean annual motion 17325693", the 
annual motion of the apse is 

148202" = 41" 10' 2". 



LECTURE XIII. 

THE MOTION OF THE APSE, AND THE CHANGE OF 
THE ECCENTRICITY. 

We have seen that when the eccentricity of the Moon's 
orbit is not considered we may write 

i = i[H-a,cos2(^-^)], 

if = ?ia» [1+ A, cos 2 (^ - ^)], 

, 3 . 2 -TO 1 , 3 TO* 
where aa=am'.-, .,- ; na = T, . 

Let us introduce the two new arbitraries e, tsr by writing 
H^hna? [1+ A5Cos(2 - 2m) d\ 

- = 7a-[l + a2Cos(2 - 2m)^ + ecos(^-'Br)], 

where A is a third arbitrary, which may be chosen to suit our 
convenience ; it must be unity when e = 0. 

Then 
dH dH dt 3 , ,^ . ^,. ^.(ft 

= - 2 m^w* ■= [sin (2 - 2m) ^ + 47ne cos (2 - 2m) e sin (^ - tsr)] 

= - 2 m^n^ , ^j [1 -(4a2 + A2)cos(2 - 2m) ^- 4e cos(^ - -or)] 
X [sin (2 - 2m) + ^me cos (2 - 2m) ^ sin (^ - isr)], 
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and also 

^^ = ^na»[l-f-A2Cos(2-.2m)fl]-na2A(2-2m)A2sin(2-2m)^. 

Now we may put A = 1 4- 17, where r) vaDishes with e. Neglecting 
powers of e above the first 

^ = ^=3m»(l + m)esin(r^2^^ + ^) 
au au 

+ 3m2(l-m)esin(3-2m^-isr)-9mS7sin(2-2m)^. 
Neglect at first the last term : 

ii = -3m«^^ ecos(r^^2m0+w)-Sm^}'~^ ecos(3^m^-isr). 
' l-2m ^ 3-2m 

Substitute this in the last term, and we get 

^ = 3m2(l+m)esin(l-.2me+tsr)+3?M«(l-m)e8in(3-2me-'Br) 



d0 



27 , 2 + 4m-4m2 . .^ . 



Now consider the other equation 

= ^'-^4-|m2nV + |m«?iV[cos(2-2m)^ 



— 2me cos (1 — 2m^ + w) + 2me cos (3 — 2m0 — tsr)]. 

Differentiate the assumed expression for -, and let h be 

chosen so that the first differential coefficient shall have the 
same form as if A, e, -bt were constant. 

Thus 
1 dr 1 /Q rt V • /r» ^ \ /\dd 1 • //I \dd 

where 

- 1 ^ [1 + a, cos(2 - 2m)^ + ^ co8(^-Br) + c ^ sin(^-i!r)= 0, 
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or 
I cos (^ - «^) + e ,^ sin (^ - -or) = 6m^ (1 + m) e sin (1 - 2m0 4- tsr) 



d0 



de 



+ 6m«(l -m)e sin (3 - 2m0-'Gr) 
2 + 4m — ^m^ 



and 



27m* 



(1- 2m) (3 -2m) 



— Grn'Oaf ^ sin (0 - isr), 



^ = 1^(2 - 2m) aasin(2 - 2m)^ + ^^ esin (0-^) 



na 



= -T- (2 - 2m)a2 8in(2 - 2m) 



na 



+ y[l+AaCos(2-2m)^]esin(^-'Br), 

g = f(2-2m)^a.cos(2-2m)^f 

-^(2-2m)A,sin(2-2OT)^esm(^-«r)^ 

+ ^ [1 + A,cos(2 - 2m) ^ eco8(^-«r)^ 
nadh,„ „ , • /a o \zi<^^ 



■'^[l+h,cos(2-2m)e][^^8m{d-w)-e'~coB(d-^)\^. 



Multiply by r*/n*a* ; then since 



r^^ = if = fena« [1+^2 008(2- 2m) ^], 



^ = (2-2m)»a2COs(2-2m)* 



we have 
nW dP 



- (2 - 2m) Aj sin ( 2 - 2m) ^ e sin (^ - tsr) 

+ e cos (0-"sr)-\- 2A2 cos (2 — 2m) e cos (0 - isr) 

-i^(2-2m)aa8in(2-2m)^ 

+ [1 + A^cos (2 - 2m)0f [^ 8in(^-'Br)-e ^ C08(^-'Br) 1 , 
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and this is equal to 
H^ u 1 r» 3 r^ . 



— 2me cos (1 — 2md + «r) + 2me cos (3 — 2m6 - -cr)] 
= H-Y + 2A3Cos(2-2m)^|[lH-aaCos(2-2m)^+ecos(^+tsT)] 

- ^3 + I m^A* [1 - 3a2 cos (2 - 2m)0 - 3e cos {6 - -cr)] 
+ 2 m«A« cos (2 - 2m) d + Qa^e cos (d - vr) 

- ^2 + 2m) e cos(l-2m^ + «•)- (g " 2^) « cos(3-2md- isr) . 

The terms in these two expressions which are independent of e 
give no new information ; equating the others : — 

- (2 - 2m) h^ sin (2 - 2m) (9 sin ((? - isr) 

-i^(2-2m)a3sin(2-2m)^ 

4-[l-f A,cos(2-2m)^pr^sin(^-isr)-e^'Jcos(tf-isr)"| 

3 3 r 

= 3m^ — ^ 7n?e cos (^ — -cr) + ^ m^ 6036 cos {0 - -cr) 

- ^2 + 2m) ^ c^s (i ^2m^ + ^) " (2 " ^^l ^ cos(3^m^-i!r) 
+ gm^iy cos (2 - 2m) 0. 
Reducing this expression 

— -^ sin (d — -cj) + e -7^ cos (^ - -cj) = ( ^ m' — - m* j 6 cos {0 - 'or) 

H-f3m«-3m' + -^m*)ecos(3-2m^-i!r) 
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and from before 
^^ //I \ 



J J 

COS (^ - «r) + e ;^ sin (^ - tsr) = 18m*e sin (d - -cr) 



+ (6m» + 6m') e sin (i - 2m^ + isr) 



+ (6m2 - 6m») e sin (3 - 2me + tsr). 



Hence, 



d5= ("4^ +-jg-m*j.sm2(^-..) 

+ ( ^m»-3m«- 3 m*^ esin(2- 2m)^ 
+ (-^^'-|^'-j|^*)^sin(2md-2isr) 

(O Q 1 P> \ 

^ m^ — g m» - Y^ mM e sin (4 - 2md — 2cr), 

dtsT 3_, 141 , , / 3 , 147 A . ^,. , 

d^ = 4^+T6"^+l 4^^ -^^r^^J cos 2(^-1.) 

+ ( ^m2-6m»+ ^m*^cos(2-2m)^ 

+ ( ^m= + |m»+ ^m^jco8(2m^-2tsr) 

+ ^-|7^a + |iyi8+ 15 m*)co8(4^^2^d-2i!r). 

We notice that among these terms, one is of long period, 
approximately semiannual, and will become of greater relative 
importance than the others on integration. 

To eflfect this integration, assume 

n = isr + asin 2 (^ - 13-)+ ^ sin (2-2m) ^ + ysin {(4 --2m)^-2i!r}, 

so that the mean motion of 11 is the same as that of vr, and 
substitute in the equation. 

A. L. 5 
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Then 

X cos 2(6 — nr) 

[Q 39 3 T 

|m«- 6m» + -g-TO* +(2-2m)^-6m»a- jm«y 

X cos (2 - 2m) 
+ r ^ "i*+ 1 »»' + fl ™' -| »»»a1 cos (2m0 - 2n) 
r 3 3 15 9 3 1 

xcos{(4-2m)d-2wj, 



so that if we take 



3 s ^219 , 

" = -8"^ +-32-'"' 

„ 9 ,^15 , 99 , 

3 . 3 , 51 , 



we have 

dn 3 

d6 



If we write mO —Il=yft, 

this becomes -^ = a — 6 cos 2-^, 



where 



3 , 309 , , 15 ,9 3^45 . 

and the solution is 

tan-i (/v/"^ b *^^ '^) = ^^^^ *' + constant. 
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Hence if we denote by -rw the mean rate of change of cj, we have 

3 , 225 , 4071 , 

d«r„ 3 , , 225 , . 4071 , 

Te= 4"^+T2"^+128"^- 

We observe that a + 6 and a — 6 are the rates of separation of 
the Sun from the apse when the Sun and the apse are at 
quadratures and syzygies with one another, respectively, — that 
is if we take n for the longitude of the apse, or, what is the 
same thing, if we ignore small terms of short period. Hence 
the mean rate of separation of the Sun from the apse is a mean 
proportional between its rates when at quadratures and syzygies 
respectively with the apse*. 

[* This is the analogue for the case of the apse of Machin and Pemberton's 
theorem on the motion of the node, inserted in the third edition of the 
Principia as a scholium to prop, xxxiii., lib. iii. See some notes by Adams 
in Brewster's Life of Newton, Appendix xxx.] 



5—2 



LECTURE XIV. 

THE LATITUDE AND THE MOTION OF THE NODE. 

Let us first treat this problem on the supposition that the 
latitude is so small that its square may be neglected. The 
equation of motion, taken from Lecture II, may be written 

where £? = r sin (latitude) and the cube of 5 is omitted; or 
neglecting the parallactic terms 






The value of /i/r* may be considered known by the operations 
which have determined the motion in an orbit coinciding with 
the ecliptic ; that is to say, 

^ = ^[1 +1^2^ H-Soa cos 2^+ ^|a22+3a4) cos 4A|rJ, 

where a has the definition of Lectures IV, V ; or numerically, 
taking 

n — n'^l, 



^ = 117 150,3 + 02523,0 cos 2t + 00025,15 cos U. 



And 

m 



^ = n'»= 00663,6. 
Hence 

^ = - ^ [117803,9 + 02523,0 cos 2< + -00025,15 cos U]. 
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Let US now consider the equation 

where P = g© + 2gi cos 2t + 2^2 cos 4it, 

in which ^i, q^ are supposed small. 

Suppose a term in z to be c cos (A;^ + /8) ; when this is 
substituted in Pz there will arise terms 



c cos (A; - 2^ + /8) c cos (A; + 2^ + )8) 



c cos (k-U + fi) c cos (A + 4e -I- /8). 

Let us therefore assume 



z = c[cos{kt-\'l3) + Ci cos{k + 2t + j3) + C2 cos(A;-h4^ + /3) + ... 



+ c_iCos(A7-2^ + )8) + c_3Cos(A7-4^ + /S)+...] 
c is arbitrary ; we have to determine k, Ci, C-i, &c. 

Substitute and equate coefficients: 

. . . + [-(A?- 4)2+9o]c-2+giC-i+9'2+ ... =0 

+giC_2+[-(A;-2)«+go]c_i+3i+?2Ci + ... =0 

+92C-2 +gic-i +[-A^+go]+giCi+g2C2+... =0 

+32C-1 + ?! +[-(A;+2)24-go]ci-+griC2+... =0 

+ 52 + qi<h +[-(A;+4)»+?o]c2+... = 0. 

If 9i> ?2>... are neglected, we have simply 

- A;^ + ?o = ; 
this is a first approximation to the value of k 

Taking qi into account and neglecting ja 

c --— 2l— 
'^ qo-{fc^2r 

c = 21. 
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In the actual case considered we notice that q^ does not 
diflfer widely from unity. Hence k is nearly equal to unity 
also, and the denominator in C-i is small, and makes c_i much 
more important than Ci. 

If we substitute these values in the third equation above, 
we have 

whence 

which may be put under the form 

whence 

With this equation we can approximate very rapidly to the 
value of k. Taking as a first approximation 

substitute this value of k in the small terms and we get as a 
second approximation 

&= 1-08516,9. 
Whence the ratio of the retrograde motion of the node to the 
Moon's mean motion is 

A/w- 1=^-1 = -00399,7, 

where g is written for A/n. This value is very correct. Taking 
the Moon's mean annual motion as 17325593", the resulting 
annual retrograde motion of the node is 

69252" = 19° 14' 12". 
Next find the values of the coefficients C-i, Ci, c_2, Cg. We have 

q^ = -01261,5, 
go -(A; -2)2= '34112,3, 
go -(A + 2)2 = -8-34022,8; 
whence as a first approximation 

c-i = - 03698,19, Ci = -00151,26. 
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Hence 

qiC-i + ga == - 00034,02, jid + 52 = -00014,49, 

and 

5o - (A - 4)« = - 7-31821, qo-{k + 4)» = - 24-6809, 

whence 

c_a = - 00004,650, 

C2= -00000,587. 

A second approximation to C-i, Ci gives 

[- (A; - 2)3 + Jo] c_i = -(?! + gaCi + giC_a), 

[-(A; + 2)« + go]ci =-(?i + g3C_i + yA); 

with the above values 

?i + q^i + ?iC-2 = -01261,47, q^ + jac-i 4- giCa = 01261,03, 

so that c_i = - -03698,00, 

Ci = -00151,20. 



LECTURE XV. 

MOTION IN AN ORBIT OF ANY INCLINATION. 

Let us consider the change in the plane of the orbit 
produced in an indefinitely small time dt by the action of a 
given disturbing force. Let Z be the resolved part of the 
disturbing force at any time in a direction perpendicular to 
the plane in which the body is moving at the instant. Imagine 
the force Z to act by impulses at the small intervals of time dt, 
then Zdt will be the indefinitely small velocity generated by 
the force Z in the time dt, in the direction perpendicular to the 
plane of orbit at the instant. 

Let FP be the radius vector and P the position of the 
body at the instant. Also let PT represent the velocity at 




the instant in magnitude and direction ; then if 2Y be taken 
perpendicular to the plane FPT and equal to Zdt, the velocity 
and its direction after the impulse will be represented by Pt, 
and the new plane of the orbit by FPt Draw Tm perpen- 
dicular to FP and join tm] then tmT is the angle through 
which the plane of the orbit has been turned about the radius 
vector FP in the indefinitely short time dt 
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Now 



tT 



tmT=:^ = ^=^ 



^Zdt 
Tm~ V 

where v is the resolved part of the velocity at P perpendicular 
to the radius vector. But 

H = vr; 

hence the angle through which the orbit is turned in an 
indefinitely short time dt is 

rZ 



H 



dt 



To find the corresponding changes in the elements that deter- 
mine the plane of the orbit, namely, the inclination of the orbit 
to a fixed plane, and the longitude of the node on that plane. 
Let NPQ be the great circle which represents the plane of the 




orbit at the time t, NR the plane of reference, usually the 
plane of the ecliptic, P the position of the body at the same 
time, and let i = PNR, the inclination, and let N be the 
longitude of the node. 

Let nPq be the position of the orbit at time t + dt 

Take NQ = 90° ; draw nm perpendicular to NPQ and qQR 
perpendicular to NR; then QR = i; and by what we have 
just proved 



NPn = ^dt 



Therefore 

where 
But 



Zr 

nm = -rfdt, sin 6, 



Zr 

qQ=jj-dt. cos 0, 



= nP. 

nm = Nn sin i = sin i dN; qQ = di 
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Therefore 

di _ Zr cos d 
dt'^ H ' 
dN ^ Zr sin 
dt H sin i ' 
which give the changes of the elements required. 

Now let NM8 be a spherical triangle, the centre of the 
sphere being 0, the centre of gravity of the Earth and Moon ; 

M 




and let 08, OM, ON point respectively to the Sun, the Moon, 
and the node of the Moon's orbit upon the ecliptic, so that NM 
is the plane of the Moon's orbit and N8 the ecliptic. Let 
M8=<o, JVS = ^, NM= 0y of which the first is identical with 
the quantity denoted by the same sjrmbol in Lecture II, but 
the second and third are not so. 

Then, following Lecture II, the forces on the Moon are 

— jr 3 COS G), in 80, 

if we ignore the parallactic terms. 
This latter may be resolved into 

^ 3 COS 0) X (cos B cos 0' + sin B sin ff cos %) in MO, 

, perpendicular to 

— ^ 3 cos G) X (sin ^ cos ^ — cos d sin ^ cos i) MO in the plane 

of the orbit, 
, perpendicular to 

— ^ 3 cos 0) X sin ^ sin i the plane of the 

orbit. 
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Now 

cos 0) = cos ^ COS ^ + sin ^ sin ff cos i 

= cos {e - ff) cos«^ + cos {0 + ff) sin^l , 
sin 6 cos ff - cos sin ^ cos i 

= sin(^- ^)cos«| + sin (^ + ^') sin^ s • 
Hence we have the following expressions for the three forces : 
^=2 + $^-i$^[{l+co«2(^-^)lcos^i + {cos2tf+cos2^} 

x2cos«|sin«|+{l4-cos2(^-d')}sin*|"|, 
r= |$[ sin2(^-^)cos^| 

+sin25.2cos«|sin2|+sin2(^+^)8in*|1, 

Z^ |^fsini[-sin(d-25')cos^| 

+ sin ^ cos i + sin (^ + 2^) sin- 1 . 
Now we have seen 

di__ Zr COS ^^_ ^rsin^ 

also, the rate of advance of the node along the orbit is 

Zr sin 
H tan i ' 

Thus the equations of motion become 

^^d^_^sin_^ 
r^~dt Hta,ni' 

together with _ = « p^ 



LECTURE XVI. 

MOTION IN AN ORBIT OF ANY INCLINATION {continued). 

To satisfy the equation at the end of Lecture XV, assume 

r = a[l+^iCos2(^-^)+^2Cos2^ + ^sCos2^+^4C08 2(d+^)], 

neglecting the square of the disturbing force and the eccen- 
tricity; thus in the small terms we write 



dO , , dff , 

r = a, -^ = n, r=a, -^ = n. 



Hence 



-^ = 7i^a[(2-2m)M,cos2(^-^') + 4^2Cos2^ 

+ imMs cos 2ff + (2 + 2m)M4 cos 2 (5 + ff)] ; 
substitute in the equation 

therefore 

H^=^ /ia[l+^iCos2(^-^) + ilsCos2^ + ^3Cos25' 

+ ^cos2(^ + ^)] 

-|mVa*r{l + cos2(^-^)}cos*| 

+ {cos2^+cos2^}2cos2|sin«^+{l+cos2(^+^)!sin*|| 

- 7i«a* [(2 - 2mf A^ cos 2 (d - ff) + 4^2 cos 2d 

+ 4m«^3 cos 2ff + (2 + 2/?i)2 A^ cos 2 (d + 5')]- 
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Again, we have the equation 
which may be written 

+ sin2d.2sin«|cos2| + sin2(^ + d08in*|l. 

Substitute for fjua its approximate value M^a* in the small terms ; 
and we find from these two equations 

- (1 - m) ^1 4- 4 (1 - my ^i 4- ^ m^ (1 — m) cos* ^ = — ^ ^' cos* ^ , 

-^2+ 4^ + 2m«.2cos^^sin«^ 

3 i . i 

= — ^ m^ . 2 cos^ ^ sin^^ , 

3 i i 

- mAi + 4m'-43 + ^ m' . 2 cos^ ^ sin'*^ = 0, 

-(l+m)^4 + 4(l+m)M44-2^Hl + m)sin*2 =--m2sin*2. 

Therefore 

. 3 , ,i 2-m 

^i = - ^m^cos*- 



2"*' ^"^ 2(l-m)(l-2m)(3-2m)' 

4 A 

il2 = — 2m2cos^^8in2^, 
ilgs Sm^cos^^sin^^ 



2 2(l-2m)(l + 2m)' 

J 3 , . ,i 2 + m 

-d.4 = — 3 m' sm* - 



and 



2"' — 2(l+m)(l+2m)(3 + 2m)' 

+ ^ :j cos* ^ cos 2 (tf - ^) H- 3m cos^ ^ sin^- cos 2^ 

3 m^ . .i «//! . /vxl 
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where a is defined by 



or 



If we preferred to define a so that the constant term in H^ 
were equal to n^a\ we should have 

n^a^ ^ fia-\- mWa^ — ^ wWa^ [cos* ^ + sin* ^ j , 

/A = nW 1 + 2 m^ - ^m^ sin^ | cos^ ^ • 

Let us next find the latitude and the motion of the node. 
Suppose that i = t'o + Ai, 

in which Ai, AJV are small, to is a constant, and Nq varies 
slowly in proportion to the time, so that we may assume 

Ai\r = i\r,8in2(d-d') + ^asin2^ + J\^8sin2^+i\^4sin2(l9 + ^), 

Ai =7i cos 2 (^-^) + /a cos 2^ + /a cos 2^' 4-/4 cos 2 (^ + d')- 
Then remembering that 

dff dN 

an expression that must be used in the terms of chief import- 
ance, we have 

^ = ^ = -2(l-m)n/,8in2(^-5')-2»/,sin2^ 

-2rm-i^)n/,8in25'-2(l + m)n/4sin2(^ + ^, 

W " ^~di ■•" ^'W" ~ ?«+2(l-»»)ni\r,co82(5-^)+2ni\r,cos2^ 
+ 2 (m-^^\ nN,<iOB2ff'+ 2{l + m)nN,cos2(e+e') 

in which the last term will be found to be required to get the 
constant q correctly to the order m\ 
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These must be equated to 

Zr cos 6 Z r sin 6 

respectively. 
Hence 

— 5 — 2mN^ — m/j -TT^ = — J m" cos i ; 

therefore as a first approximation 

g 

hence 

0/1 \ r 3 o • • o^' r 3 m* . . ,i 

— 2(1 -m)/i= -rm^smicos*^, /i= — -, smtcos*^, 

^ ^4 2 81— m 2 

3 . 3 . . . 

— 2/2 =— ^m^sinicosi, /2= ^ m* sin i cost, 

3 „ A ^ 3 „ . . r- 3 m sin i 



(3 \ 3 

m+ -m^cosi j /g = — ^m^sin i, /8= 



1 + T ^^ cos i 

4 



-2(l+m)/4 = --m«smtsin22' ^4= ST+m®^^* 2' 
and 

3 3 

2i\ra = 7 m* cost, i\ra= o^^cosi, 

4 o 



/ 3 \ 3 3 

2fm+Tm*cosijiV8= 7m*cosi, N^= ^ 



3 mcosi 



3 * 

1 + 7 m cos i 
4 



3 m* 



Substitute above for the quantities /j, N^ and we get the second 
approximation to 5, 

3.9 9 

5 = jm* cos i — ^ m* cos* * + ^ ^' sin*i. 
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It will be observed that /g, N^ are of lower order than the other 
coeflScients, so that in order to obtain them correctly to the 
same order as the others we were obliged to retain small terms 

in -^ arising from the variability of N, 

If we take the variable plane defined by the longitude of the 
node Nq and the inclination % as the plane to which the position 
of the Moon is referred, we have the latitude of the Moon above 
this plane 

= Ai sin 6 — A^ sin i cos 

1 



= Q m sin I cos^ ^ 



m 



1.3 

1 + T ^ COS I 

4 



1 —m 



sin (5- 219') 



— 5 m* sin % cos i sin 

o 



3 . . . J 

+ 5 m sm 1 sm^ ^ 

o Z 



m 



1 + T wi cos i 
4 



l+m 



sin (5 + 219'). 



LECTURE XVIL 
ON hill's method of treating the lunar theory. 

Let us suppose the Moon to move in the plane of the 
ecliptic, and let us refer its motion to rectangular axes in 
rotation, the rotation being such that the axis of x passes 
always through the mean position of the Sun; that is, the 
axes rotate with angular velocity n\ and if we suppose the Sun 
describes a circular orbit about the origin, its coordinates are 

x'^a\ 2^ = 0. 

Let X, y be the coordinates of the Moon. 

Then the disturbing forces of the Sun upon the Moon 
relative to the Earth are 

rn! x — a' rnf m' y 

parallel to the axes of x and y respectively, where 

p'^ix^ay + j/', 

and the forces of the Earth on the Moon relative to the 
Earth are 







/JL X 


r' r' 


where 




r« = 


c^ + y\ 


Now these forces 


may 


be written 






dn 


dil 






dx 


' dy' 


where 




r 


m' m'x 
■^ 7 ~ "a'» 


A. L. 
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Hence 

We have tacitly assumed the origin to be at the centre of the 
Earth ; if we prefer to place it at the centre of gravity of the 
Earth and Moon, the necessary change is effected by multiplying 
the last terms, which correspond to the Parallactic Inequalities, 
hy{E-M)l{E+M). 

Equating these forces to the accelerations of the Moon 
parallel to the coordinate axes, we have the equations of 
motion in the form 

dt" ^"^ dt ""^ da>' 

or, as they may be written, 

dt' dt dx' 

dt" ^ di ~ dy ' 
where i2 = ft + 5 n'" («» + y) 

Now suppose we have found values of x and y which satisfy 
this pair of equations and which involve two arbitrary constants. 
This may be accomplished by taking assumed developments 

X = Soi cos i {t + 7), 

y = %hi sin i {t + 7), 
substituting in the equations, and equating coeflBcients of the 
various terms. The solution found will include the Variation 
and the Parallactic Inequalities. Let it be required to amend 
this solution by the introduction of the remaining two arbitrary 
constants that are required for a complete solution. 
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Let the additional terms that we seek be Sx, Sy, which we 
shall suppose so small that their squares and products may be 
neglected, let us consider first the terms which are multiplied 
by the first power of one of the new arbitraries, the original 
particular solution corresponding to the case in which this 
arbitrary is zero. 

Then Sx, Sy are determined by the equations 

d^Sy ^ ,dSx d^R ^ d^R ^ „ 

dt"" dt dxdy dy^ ^ 

where X, Y are supposed known functions of x, y or of t, and 

have been added here to include disturbing causes not allowed 

for in the above form of R. 

dx dii 
Multiply the original equations by -it , -^ and add : 

d?x dx d^y dy __ dR dx dR dy 
d^diWdi'~dxdi^dydi 

_dR 

~ dt ' 
since x, y are the only functions of t that R involves ; whence 



(S)'-(l)'=^^-<'. 



where G is an arbitrary constant ; this is the integral known as 
Jacobi's Integral. 

Let us write 

then we have V^ = 2R + C; 

and from the original equations themselves 

^=^cos0 + ^sm^ 

dR , . dR . - 
= ,- cos <p + T- sm <p ; 
dx dy ^ 
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and 

dx dxi 
And from these, differentiating and substituting for ^, -^, 

we get 

dt^ ^ dt U^ J~ L d^ ^ 

+ ^-d^''''"f'''^'f'+df'''''^\' 
^ d^<f> . „ dV (d^ , ,\ „V d-'R . . . 



LECTURE XVIII. 

ON hill's method of treating the lunar 
THEORY (continued). 

The equations for &c, By are 

the equations for x, y are 

cP^_„ ,dy _dR 
dP~ di~'dx' 



?l + 2n 






Multiply the former pair ^y -ji^ ~^ respectively, and the 
latter pair by -rr > ^> and add all together; we get 

dxd^ix d^x dSx dy d^Sy d^y dhy 
m 'djf^li^'df'^'M'df'^'dt^llJt 

(d^Rdx^R dy\^ f d^R dx d'R dy\ ^ 
" \da? dt ^ dxdy dtl ^^ "*" \dxdy dt "*" dy* dt) ^ 

dR dSx dR dSy y^^ , y^^ 
dx dt dy dt dt dt ' 

Now (PRdx dlB dy^ddR 

dx' dt dxdy dt dt dx* 

d^R_ dx d^dy _d dR 
dxdy dt dy^ dt " dt dy' 
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Then our equation may be integrated 

SO that r is a known function of t, which involves an arbitrary 
constant. 

Now let us assume 

Ba; = v cos ^ — w sin <^, 
Sy = V sin ^ + w cos ^. 

_.,,.. , ,. dx dy dBx d8y „ , 

Substitute above ^o^ 7s > ;£ > jT > df''^^ 

„(dv dd>\ fdR . , dR . .\ 
^[dt-''dt) = [d^'''"f'^^''''V' 

( dR . . dR ^\ , ^ 

dR ^ dR . ^ dV 
But _cos,^ + -^sm,/, = -^, 

Therefore F(f-.f) = f.+ F(f + 2«')..r, 

whence •=.= /= (-^ + n'j wdt + 1 p^ d^. 

An arbitrary constant is included on the right. This 
equation shews that when w is known, v can be found; it 
remains to determine w. 

Now by actual diflferentiation 



, dhx . . , diy dv d6 
cos,^-^ + sm0^ = ^-«;^ 



. ^ d^Sx ^ d»Sy d^w „ dv d4> /diV , 



cP<f> 
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Also multiplying the differencial equations for Sx, Sy by 
— sin (f>, cos <^, respectively and adding 



— JT sin ^ + Fcos ^. 

Substitute and we find 
dhu „dvdA (d4>\^ . d^6 , ^ , /dv dd>\ 

\ d?R . . .^ d^R . ,. . ...^d^R . . .1 
" ^[— ^sin(^cos^+ ^^(cos'^-sin»^)+-^-,8m^cos<^J 

V d?R . .. .d'R.. . .d^R..! 

— Xem<f> + Yoos<l>. 

Now we have seen 

dv vdV^^fd<l>^ A j.^ 

- ^-v^+^Kdt+^'r'-v 

Substitute for 2 f ^ + n'j -^ on the left. 
We get 

"-[S4f(t«')]-Kt-')'-©'-'t] 

[ d'R . . . , d?R, .. . , .. d'R . . r\ 

— X sin ^ + Fcos <^. 



dhv 
dt 



88 LECTURES ON THE LUNAR THEORY. [LECT. XVIII. 

But by the equations proved at the end of Lecture XVII, 
the terms in v cancel one another, and we are left with the 
equation for w: 

d'R 
dxdy* 

= -2(^^ + w')^-Z8in<^+Fcos<^. 

Or since 






. 1 da: . , I dy 
#^^ , lfdR,.dR A 

the coefficient of w is 

^(_dRdy dR d^y __a!}Lf_^^ dRdx\ ,^ 
V\ dxdt'^dydt) ^ kA dx dt'^ dydt)'^ 

F^ I dar^ U^/ da;dy de d^ "^ dy^ \dt) \ 
= P, say. 

This function P is a known function of ^ ; it may be seen 
that if 

X = 2ai cos i {t + 7), 

y = 26»8ini(]t + 7)> 
then P may be developed in the form 

P = 2il» cos i (^ + 7). 

Hence if we omit the terms X, F, due to other disturbances 
not yet allowed for, the equation for w assumes the form 

^ + w [^0 + ^icos (^+ 7) + ^2 cos 2 (^ + 7) + ...] = 0. 

This is identical in form with the equation treated in 
Lecture XIV, to find the motion of the node. The value of w 
may be found by the method there employed, and the value 
of V deduced from it. 
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